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PARTITION FUNCTIONS AND THERMODYNAMIC PROPERTIES 


TO HIGH TEMPERATURES FOR H3 AND H^ 
by R. W. Patch and Bonnie J. McBride 
Lewis Research Center 

SUMMARY 

Tables of partition functions were compiled for Hg and Hg at temperatures from 
298. 15° to 56 000° K. Tables of thermodynamic properties were compiled at temper- 
atures from 298. 15° to 10 000° K. The latter tables give the following thermodynamic 
functions for ideal gases: heat capacity at constant pressure C^/R, sensible enthalpy 
(Ht - H 0 )/RT, entropy at 1 atmosphere S^,/R, sensible free energy at 1 atmosphere 

-(G^, - Hq)/RT, enthalpy H t /RT, and free energy at 1 atmosphere -G^/RT. The heats 
of formation at 298. 15° K are also given. Since no band spectra have been observed for 
either ion, all calculations were based on ab initio potential energy calculations found in 
or derived from the literature. All results for Hg are reliable up to 15 000° K, but 
results for Hg are tentative at all temperatures. 

A digital computer program to calculate Wentzel -Kramers -Brillouin rotational 
energies from the potential energy of a diatomic molecule or molecular ion is also 
included. The program is in FORTRAN IV. 

INTRODUCTION 

In very -high-temperature propulsion devices such as gaseous nuclear rockets, the 
dominant mechanism of heat transfer is radiant energy exchange between volumes of 
plasma, and between plasma and the wall. To calculate such heat transfer, it is 
necessary to know the opacity of the plasma. For a hydrogen plasma under many condi- 
tions, minor species make major contributions to the opacity, either directly or by their 
influence on the composition. Two such minor species are the hydrogen triatomic 
molecular ion Hg and the hydrogen diatomic molecular ion Hg. The partition functions 
of these two molecular ions are needed to calculate the composition of a hydrogen plasma 
at high pressures and low degrees of ionization, and the thermodynamic properties are 



needed to interpret data from some experimental investigations of the opacity of hydrogen. 

The existence of Hg and H^ has been well established by investigations using mass 
spectrometers, but no band spectrum of either ion has been observed. There have been 
no published calculations of the partition function or thermodynamic properties of Hg, 
and only one such calculation for Hg (ref. 1) has come to the authors' attention. 

The purpose of the present work was to calculate the partition function for tempera- 
tures between 298. 15° and 56 000° K and thermodynamic properties for temperatures 
between 298. 15° and 10 000° K for Hg and Hg. This presented two principal difficulties: 

(1) complete sets of spectroscopic constants were not available because no band spectra 
have been observed, and (2) the usual approximations used in calculating partition func- 
tions are not accurate at high temperatures. 

The solution was to utilize existing ab initio calculations for the potential energy of 
both molecular ions. For Hg the potential energy was used to obtain both low- and high- 
temperature approximations to the partition function; the two approximations were then 
faired together. For Hg the rotational term values were calculated from the potential 
energy by the Wentzel -Kramers -Brillouin (WKB) approximation, and the vibrational 
term values were obtained from more accurate existing calculations. The partition 
function was then found by direct summation over the vibration-rotation states. 

The present work was limited to ideal gases with no excited electronic states. 


TRIAT0MIC MOLECULAR ION H 3 

The calculations in this section are divided into two sections, the Analysis and the 
Results and Discussion. The latter of these contains numerical results for the partition 
function and thermodynamic properties and a discussion of the accuracy of the partition 
function. 


Analysis 

Assumptions and re strictions. - The assumptions and restrictions utilized in the 
analysis are listed here for convenience. They are the following: 

(1) The Born-Oppenheimer approximation (ref. 2) is assumed to be applicable. 

(2) The equilibrium configuration is an equilateral triangle. 

(3) There are no excited electronic states (this is a good assumption, since the 
lowest excited electronic state is about 75 000 centimeters”* above the ground state and 
appears to be a repulsive state (ref. 3). 
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(4) Differences in nuclear spin degeneracy are neglected in computing the thermo- 
dynamic properties. 

(5) An ideal gas is assumed. 

Potential energy . - A number of ab initio calculations (refs. 3, 4, 5, 6, and 7) and 
one semitheoretical calculation (ref. 8) have been made of the potential energy of Hg. It 
is now generally agreed that the equilibrium configuration is an equilateral triangle with 
a proton at each vertex. The results of these calculations are given in table I. All but 
the earliest (ref. 3) ab initio calculations give about the same equilibrium internuclear 
distance r e , but there is considerable disagreement about the equilibrium potential 
energy V . (All symbols are defined in appendix A. ) 

There are no experimental values for the dissociation energy or potential energy of 
Hg or any of its isotopic modifications, but Stevenson and Schissler (ref. 9) measured 
the rate of the reaction 


D 2 +D 2- D 3 +D 

where Dg is the deuterium diatomic molecular ion, Dg is the deuterium molecule, Dg 

is the deuterium triatomic molecular ion, and D is the dueterium atom. It was concluded 

that the rate of this reaction was so high that the reaction was unlikely to have any 

activation energy (private communication from David P. Stevenson, Shell Development 

Company, Emeryville, California). An activation energy of zero gives a lower limit to 

the dissociation energy of D« (into D + D + D + where D + is the deuterium atomic ion) 

- 1 J 

of 58 456 centimeters . In the Born-Oppenheimer approximation the potential energies 

of Hg and Dg are the same. If a Dg zero -point energy of 3684 centimeter - * calculated 

from Conroy’s (ref. 5) Hg potential energy is assumed, the equilibrium potential energy 

of Dg (and hence Hg) is thus -6. 214x10^ centimeter - * or less. This value is consistent 

with all calculations for Hg except Hoyland's (ref. 6). 

Conroy's equilibrium potential energy (ref. 5) is generally considered to be the most 

reliable for three reasons: (1) Conroy used a 26-term wave function, (2) Christoff ersen 

(ref. 4) estimated that if he included more terms in his calculation, he would get an 

4 -1 

equilibrium potential energy of -7.68x10 centimeter , which is close to Conroy's value, 
and (3) Conroy and Bruner (ref. 10) calculated the potential energy of the linear hydrogen 

triatomic molecule H„ by a development of the method used for Ho and found the rim of 

« 3 -1 J 

the potential energy basin of Hg to be 2. 17x10 centimeter above the energy of Hg + H 

where Hg is the hydrogen diatomic molecule and H is the hydrogen atom. This is in good 

agreement with the experimental activation energy (ref. 11) for the exchange reaction 

H + Hg — Hg + H 
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Consequently, Conroy’s Hg potential energy was used in this report. 

As explained in appendix B, it was necessary to make a least-squares fit to Conroy’s 

Hg potential for use in the low -temperature approximation. This resulted in slightly 

different values of r and V which are also given in table I and which for consistency 

6 0 

were used in calculating the partition function and thermodynamic properties in the low- 
temperature approximation. 

Partition function. - The energy of the rotationless ground vibrational state was used 

as the reference for the internal partition function qTT+, which was obtained by two dif- 

n 3 

ferent approximations. For low temperatures it was computed from the calculated 
spectroscopic constants (eq. (B39)) by means of Woolley's approximation (ref. 12). For 
high temperatures a classical partition function with vibrational quantum corrections 
(eq. (C21)) was used (see appendix C). 

Thermodynamic pr operties. - The thermodynamic functions C p /R, (H t - Hq)/RT, 
-(G^, - Hq)/RT, and S^,/R for ideal gases were calculated on a digital computer. Equa- 
tions for evaluating the thermodynamic functions from the logarithm of the partition 
function and its derivatives are as follows (ref. 13, ch. V): 

C p _ t 2 d 2 (ln q) ( 2T d(ln q) [ 5 ^ 

R dT 2 dT 2 

H T ~ H 0 = T d(ln q) + 5 ( 2 ) 

RT dT 2 

s° 

_T = T d(ln q) + ln q + 1 i n M + 1 In T - 1. 164953 (3) 

R dT 2 2 

-(G° - H n ) S° H t - H n o 5 

T 2_ = _£-_£ 2 = In q + .“In M + 2 In T - 3.664953 (4) 

RT R RT 2 2 

where M is the molecular weight or 3. 02336 grams per mole for Hg. Physical constants 
have been taken from reference 14. 


Results and Discussion 

Partition function. - The low-temperature partition-function approximation with 
first-order corrections (ref. 12) and with both first- and second-order corrections are 
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shown in figure 1. Above 7000° K the second-order correction was appreciable compared 
to the first-order correction, causing one to suspect that higher-order corrections 
should be included. In addition, at high temperatures the least-squares fit (eqs. (B15) 
and (B17)) used to calculate the spectroscopic constants is inadequate: it does not fit 
Conroy’s potential energy at the high energies corresponding to some of the energy levels 
which contribute significantly to the partition function. 

The high-temperature approximation to the partition function q„+ is also shown in 

figure 1. For low temperatures, it gave about twice the values that the low -temperature 
approximation did. The principal reason for this is the crude but necessary assumption 
in the high -temperature approximation that the potential energy V is the sum of three 
two-atom potential energies U (eq. (C 5)) coupled with the assumption that U is one- 
third of Conroy’s (ref. 5) V for the equilateral triangle configuration. This results in 
constant-V lines that tend to fall outside Conroy's constant-V lines (fig. 2) for isosceles 
triangles, except, of course, along the line Q corresponding to equilateral triangles. 
Similar discrepancies would undoubtedly occur for many scalene triangles if Conroy had 
calculated V for scalene triangles. Hence, the high -temperature approximation to the 
partition function gave values that were too large at low temperatures. 

At temperatures of the order of 40 000° K the high -temperature approximation con- 
tains at least two significant inaccuracies in addition to the one mentioned in the 
preceding paragraph: (1) all metastable rotational states are omitted and (2) the 
approximation yields the following identical heats of reaction: 

Hj - H+ + H AH q = 52 200 cm -1 

Hg H 2 + H + AH q = 52 200 cm -1 

(using potential energy for Hg from ref. 5) whereas the best available information 
(refs. 5, 15, 16, and 17) used directly gives 

Hg — Hg + H AH q = 55 800 cm" 1 

Hg - H 2 + H + AHq = 40 100 cm" 1 

(zero-point energies have been neglected in all four equations). Inaccuracy (1) tends to 

make q H + too low whereas inaccuracy (2) and the inaccuracy in the preceding paragraph 
n 3 

tend to make qxr+ too high. Thus at temperatures of the order of 40 000° K the inac- 
curacies in the high -temperature approximation tend to cancel. 
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Taking the inaccuracies mentioned in the preceding three paragraphs into consider- 
ation, a smooth curve (fig. 1) was used to fair the low-temperature approximation into 
the high -temperature approximation. This curve extended from 7000° to 26 000° K and 
had the formula 

In q H + = 3. 1705 + 1. 03824xl0" 3 T - 3. 2449xl0~ 8 T 2 

+ 5. 2127X10 -13 T 3 - 4. 077xl0" 18 T 4 (5) 

The resulting q^+ is shown in figure 3 and is tabulated in table n. 

Thermodynamic pro perties. - The thermodynamic functions C p /R, (H T - Hq)/RT, 
S°/R, -(G° - H Q )/RT, H t /RT, and -G°/RT for Hg are tabulated in table IH. The 
functions are given for temperatures ranging from 298. 15° to 10 000° K at 100° K 
intervals from 300° K. The absolute values, H t /RT and -G^/RT, are relative to an 
assigned enthalpy base of Hggg ^ = 0 for Hg and e. 

The following heats of reaction AHq were used to obtain the heat of formation of Hg 
from Hg(g) and electron gas e: 

Hg - 2H + H + AH q - 72 707 cm" 1 (appendix B) 

Hg — 2H AHq = 36 113.6 cm" 1 (refs. 16 and 18) 

H — H + + e AHq = 109 678. 76 cm" 1 (ref. 19) 

where H + is the hydrogen atomic ion. At 298. 15° K the following values of (Hggg 
- Hq)/RT were used: 3. 41561 for Hg (ref. 20), 2. 5 for e, and 3. 9745 for Hg (tabie HI). 

The heat of formation at 298. 15° K was calculated to be: 


| H 2 - H 3 + e 


AHggg 15 = 1 093 640 J/mole 
= 261 386 cal/ mole 


Accuracy . - The estimated accuracy of q H + (table H) is different at different tem- 
peratures. Assuming Conroy’s V is exact, qrr+ should be accurate to within 20 percent 

n 3 

from 298. 15° to 8000° K. From 8000° to 15 000° K, qTr+ should be accurate to within 

H 3 

a factor of 2. Above 15 000° K, qrr+ is an order of magnitude estimate. However, at 

n 3 

1000-atmospheres pressure and above about 8000° K there is less Hg than Hg or H + in a 
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hydrogen plasma. (Below about 8000° K there is more Hg than Hg and at still lower 
temperatures there is more Hg than Hg or H + ). At lower pressures, the temperature 
where the concentrations of H+ and H + are equal will be less than 8000° K. Hence, 
errors in q + at high temperatures are not of much practical importance. Considering 

the almost complete lack of experimental data on Hg, a more sophisticated calculation 
of qjj+ than is contained in this report is not warranted. 


DIATOMIC MOLECULAR ION 

The calculations in this section are divided into two parts, the Analysis and the 
Results and Discussion. The latter of these contains numerical results for the rotational 
term values, partition function, and thermodynamic properties, comparison with other 
investigators, and a discussion of the accuracy of the results for Hg. 


Analysis 

Assumptions and restrictions . - The assumptions and restrictions utilized in the 
analysis are listed here for convenience; they are the following: 

(1) There are no excited electronic states (this is a good assumption below 15 000° K 
because the lowest electronically excited state is repulsive and the next excited state is 
about 92 600 centimeters - '*' above the ground state and only has a shallow potential well 
(ref. 21)). 

(2) Differences in nuclear spin degeneracy were neglected in computing the thermo- 
dynamic properties. 

(3) An ideal gas was assumed. 

Potential energy and energy levels . - Experiments and various ab initio calculations 
of the constants of Hg have resulted in the values listed in table IV. The experimental 
values of r g and u>' have some uncertainty because no band spectra of Hg have ever 
been observed. Obviously, the small differences in the various values for the constants 
in table IV would have negligible effect on the partition function and thermodynamic 
properties, and it can be safely assumed that the potential energies used by Cohen, et al. 
(ref. 22) and by Wind (ref. 15) are essentially correct. 

The best potential energy that could be obtained from the literature was found by 
adding Wind’s electron energy (ref. 23) to the nuclear repulsion potential energy. Rota- 
tional term values were obtained from this potential energy by the WKB connection 
formula for Hund's coupling case (b) (eq. (D10)) with the help of equations (Dll), (D12), 
and (D16). The second term on the right hand side of equation (D12) was assumed to be 
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zero. The computer program used for this calculation is described in appendix D. 
Vibration -rotation term values &\v, K) were obtained by adding vibration term values 
from Wind (ref. 15) to the rotational term values obtained from the computer program. 
Partition function. - The energy of the rotationless ground vibrational state was used 

which was obtained from the 


as the reference for the internal partition function q 
relation 


H + ’ 

n 2 


q H , = ^ Y Y (2K+ l)e- hc ^- K >/ kT <«> 

v K 

where g g is the electronic statistical weight, a is the symmetry number, v is the 
vibrational quantum number, K is the quantum number for total angular momentum 
apart from spin, h is Planck's constant, c is the velocity of light, k is Boltzmann's 
constant, and T is temperature. Metastable rotational levels were included in the 
summation in equation (6). The quantities g g and a are both equal to 2. In equation (6) 
the difference in statistical weights of ortho and para states is neglected as was done for 
Hg. This is a good approximation at 298. 15° K and above. 

Thermodynamic properties. - The thermodynamic functions C^/R, (Hrp - Hq)/RT, 
-(G^, - Hq)/RT, and S°/R for were calculated on a digital computer using equa- 
tions (1) to (4). The molecular weight was taken to be 2. 01539 grams per mole. 


Results and Discussion 

Rotational term values. - The calculated rotational term values are given in table V. 
Note that these rotational term values are not the same as obtained in the sample problem 
(tables VI and VII) because many more potential energy points with internuclear distance 

Q 

r as great as 25x10 centimeter were provided as input to the program when table V 
was calculated. Where no value is given in table V, the state is nonexistent. The rota- 
tional term values in table V agree within 0. 85 centimeter”* with rotational term values 
calculated from Wind's energy levels (ref. 15) and within 0. 12 centimeter"* with term 
values from reference 24. 

Partition function. - The partition function calculated from equation (6) is given in 
table II and figure 3. It agrees within 38 percent with the less accurate values of Vardya 
(labeled r^^ (nj = 40) in table III of ref. 1), who made calculations for 3000° to 
25 000° K. 

Thermodynamic properties. - The thermodynamic functions Cp/R, (H^, - Hq)/RT, 
S°/R, -(G^, - Hq)/RT, H^/RT, and -G°/RT for Hg are tabulated in table VUI. The 
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functions are given for temperatures ranging from 298. 15° to 10 000° K at 100°K 
intervals from 300° K. The absolute values, H t /RT and -G^/RT, are relative to an 
assigned enthalpy base of Hggg 15 = 0 for Hg and e. 

The following heats of reaction AHq were used to obtain the heat of formation of Hg 
from Hg(g) and electron gas e: 


H 2 -2H 

AH Q = 36 113.6 cm" 1 

(refs. 16 and 18) 

H — H + + e 

AH q = 109 678.76 cm' 1 

(ref. 19) 

Hg — H + H + 

AHq = 21 379.36 cm" 1 

(ref. 15) 

Using (H 298 15 - H Q )/RT values of 3. 41561 for Hg (ref. 20), 2. 5 for 
Hg (table VUI) results in the heat of formation at 298. 15° K becoming 

e, and 3. 4623 for 

H 2~ H 2 + e 

AH 298 15 = 1 610 J / mole 



= 357 220 cal/mole 



Accuracy. - The estimated accuracy of qrr+ (table n) is different at different tem- 

2 

peratures. The contribution of the ground electronic state is given to within 1 percent by 

equation (6) for 298. 15° to 56 000° K. However, at temperatures above 15 000° K there 
are contributions of 1/ 2 percent or more due to electronically excited states with 
potential wells (ref. 21) which are not included in equation (6). At 56 000° K the value of 
q H + in table n is consequently about 70 percent too low. The large error at high temper- 
atures is of little practical importance since at 1000-atmospheres pressure and 15 000° K 
there is only about 4 percent as much Hg as H + in a hydrogen plasma and even less at 
higher temperature or lower pressure. 


CONCLUDING REMARKS 

In this report the partition function and thermodynamic properties of Hg and Hg were 
calculated for ideal gases. Because no band spectra have been observed for Hg or Hg, 
potential energies had to be obtained from ab initio calculations found in the literature. 

The partition function and thermodynamic properties of the ground electronic state 
of Hg are based on reliable potential energies and other calculations and are unlikely to 
change appreciably if the calculations are carried to higher degrees of approximation or 
if more experimental information becomes available. However, if excited electronic 
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states of Hg were included, the partition function for would be appreciably higher at 
temperatures in excess of 15 000° K. 

The partition function and thermodynamic properties of Hg are based on the most 
reliable potential energies available from ab initio calculations, but there is very little 
published experimental information. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 18, 1967, 

122-28-02-17-22. 
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APPENDIX A 


SYMBOLS 


A i 


B . C 
e’ e 

Bj(i=l, 2, . • , 7) 

B , C 
v’ v 

c p 

c 

D J’ D Jk ,D k’ p 

D o 

E 

E’ 

E F 

E o 

Ffry) 

•*(v,j) 

^ t (v, j) 
f(v,j) 

G 

G -1 

Grp 

Se 

H t 


nondegenerate totally symmetric symmetry species for point group 
Dg h (ref. 25) 

h 2 /8jr 2 pr 2 

rotational constants for equilibrium separation of atoms 

constants in expression for potential energy in terms of £ 

rotational constants for vibrational state v^, Vg 

heat capacity at constant pressure 

velocity of light 

rotational stretching constants 

energy to dissociate rotationless ground vibrational state 
energy 

doubly degenerate symmetry species for point group (ref. 25) 
vibration-rotation energy eigenvalue 
unevaluated constant term in energy 

ratio of incomplete gamma function to complete gamma function defined 
in ref. 26 

vibration-rotation t erm value for actual potential (an eigenvalue if 
corresponding v has the value 0 or is a positive integer) 

vibration-rotation term value for Morse potential 

rotational term value 

kinetic energy matrix (ref. 25) 

inverse of kinetic energy matrix (ref. 25) 

Gibbs free energy at temperature T 

electronic degeneracy 

sum of sensible enthalpy at temperature T and chemical energy at 
0° K 


H t - Hq sensible enthalpy 
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AH t 

H o 

3C 


j(e) 

ii 


Ii 


j 

K 


k 

k lll’ k 122’ k 

k 1122’ k 2222 

L 


1111 ’ 


/ 


M 

m 

P(r) 

Pl’P2’P3 

Q 

Q.(i=4, 5, . . . , 9) 

Q 1 ’ Q 21’ Q 22 

q 

q c 


heat of reaction at temperature T 
chemical energy at 0° K 

Hamiltonian of Hg with reference energy of infinitely separated 
H + H + H + in ground states 

Planck’s constant 

equilibrium moment of inertia about principal axis i 
integral due to two -bond states of Hg in expression for Z’ 
integral due to three -bond states of H^ in expression for Z’ 
total angular momentum quantum number 
generalized rotation quantum number 

quantum number for total angular momentum apart from electron 
spin 

Boltzmann constant 
force constants 

electronic orbital angular momentum vector 

matrix for orthogonal transformation from mass -weighted 
Cartesian coordinates to normal coordinates (ref. 25) 

vibrational angular momentum integral number 

molecular weight 

mass of proton plus two -thirds of mass of electron 

radial wave function of diatomic molecule or molecular ion 

momentum of atoms 1, 2, and 3, respectively 

normal coordinate column vector for Hg 

normal coordinates for rotation and translation of 

normal coordinates for Hg vibration 

internal partition function 

classical internal partition function of H^ including only bound 
states and with reference energy of rotationless ground vibra- 
tional state of Hg 


12 



IH? I®« 


<l'o 


<*Ct 


^Ct 

R 

R,R' 


r ’ r ij’ r 12’ r 13’ r 23 

r e 

r max 


r min 


r 0l’ r 02 

r 03 

r 04 

r v r 2 ,r 3 


S 

S T 


pg 2’ <^3 


S 1 ’ S 21’ S 22 

T 


U ’ U a’ U b’ U o 

U a’ U b’ U a 


classical internal partition function of Hg including only bound 
states and with reference energy of infinitely separated 
H + H + H + in ground states 

classical partition function of Hg including translation of the mol- 
ecular ion and only bound states and with reference energy of 
infinitely separated H + H + H + in ground states 

classical partition function of Hg including translation of the atoms, 
bound and unbound states, and with reference energy of 
infinitely separated H + H + H + in ground states 

universal gas constant 

internuclear distances (fig. 2 only) 

ratios of quantum -mechanical harmonic -oscillator partition func- 
tion to classical harmonic -oscillator partition function 

internuclear distances 

equilibrium internuclear distance 

internuclear distance of far turning point 

internuclear distance of near turning point 

internuclear distances of inflection points in U” (fig. 6) 

internuclear distance of minimum of U M (fig. 6) 

internuclear distance of maximum of U" (fig. 6) 

position vectors of atoms 1, 2, and 3, respectively, relative to 
fixed axes 

total electron spin quantum number 
entropy at temperature T 
£ column vector 
£ row vector 

changes in internuclear distances of Hjj 

symmetry coordinates of H^ 
absolute temperature 

diatomic potential such as obtained from RKR approximation 
diatomic electronic potential (including nuclear repulsion) 
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U”,U” 

effective diatomic potential 

u t 

diatomic Morse potential 

U 12’ U 13’ U 23 

two -atom potential energies of Hg with reference energy of 
infinitely separated H + H + IT 1 " in ground states 

U Z 

dimensionless energy of fundamental vibration l of Hg 

V 

three -atom potential energy of Hg with reference energy of 
infinitely separated H + H + H + in ground states 

V’ 

three-atom potential energy of Hg (zero for equilibrium inter- 
nuclear distances) 

V e 

three-atom equilibrium potential energy of Hjj with reference 
energy of infinitely separated H + H + H + in ground states 

v,v r v 2 

vibrational quantum number if integer; vibration number if 
noninteger 

V 

max 

maximum value of v 

W 

volume of container 

W a 

diatomic electronic potential energy for Hund's coupling case (a) 

X 

Cartesian coordinate 

X 7 7 
L 2 L 2 

Ax : , AX 2 ,AXg 

vibrational angular momentum constant 

displacements of protons 1, 2, and 3, respectively, in x-direction 

X 11 ,X 12 ,X 21’ X 22 

anharmonicity constants 

y 

Cartesian coordinate 

Ay r Ay 2 ,Ay 3 

displacements of protons 1, 2, and 3 respectively, in y- 


direction 


Z’ 

Z" 

z 


Hg configuration integral including only bound states 
Hg configuration integral including bound and unbound states 
Cartesian coordinate 


z' 


AZj, AZg, AZg 

a 0’ a l’ a J 

y 0’ r l ,r 2 
r(s) 


Cartesian coordinate with origin at proton 1 

displacements of protons 1, 2, and 3, respectively, in z -direction 

vibration -rotation interaction constants for Hg 
representation 
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I 


h 

v 

7j. (i=l, 2, . . . , 9) 

9 1 ,9 2 

K 

A 

X m’ X V X 2 

M 

S 


O’ 

<Pl><P2 

* 

ft 


Subscripts: 

H 2 

•*3 

j 

k 

T 

Superscripts: 


Coriolis splitting constant for Hg 
r i column vector 

mass-weighted Cartesian displacement coordinates of Hg 
polar angles of Hg (fig. 5) 

quantum number for component of total angular momentum along 
z axis for 

quantum number for component of orbital angular momentum 
along internuclear axis 

squares of angular velocities of normal vibrations of Hg 

reduced mass of two atoms in diatomic molecule or molecular ion 

quantum number for component of electron spin angular momentum 
along internuclear axis for diatomic molecule or molecular ion 

symmetry number 

aximuthal angles of Hg (fig. 5) 

vibrational wave function of diatomic molecule or molecular ion 

quantum number for component of total electronic angular momentum 
along internuclear axis for diatomic molecule or molecular ion 

zero-order vibrational frequencies of Hg 

fundamental vibrational frequencies 


hydrogen diatomic molecular ion 
hydrogen triatomic molecular ion 
proton j of Hg 
vibrational kinetic 

temperature given as table argument 


0 
* 

1 


time derivative 

standard state (1-atm pressure) 
effective value for unbound atoms 
effective value for bound atoms 
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APPENDIX B 


SPECTROSCOPIC CONSTANTS AND DISSOCIATION ENERGY OF H3 

Since the spectrum of Hg has never been observed, it is necessary to obtain the 
spectroscopic constants and dissociation energy from theoretical calculations. This is 
done in this appendix, using the potential computed by Conroy (ref. 5). 

The procedure consists of three parts: (1) The vibration of the molecular ion is 
analyzed assuming that the potential energy V contains only terms that are quadratic in 
the symmetry coordinates or normal coordinates. This is known as the harmonic oscil- 
lator approximation and gives zero-order frequencies of vibration, which are the 
frequencies which the molecule would exhibit if the anharmonic terms in V were 
all zero. (2) The symmetry of the molecular ion is used to find the form of the 
anharmonic terms in V, and a least-squares fit to the potential computed 
by Conroy is made. (3) The anharmonic terms in V are treated as perturbations 
to the harmonic oscillators. With the exception of the least-squares fit, this pro- 
cedure has been well established for molecules (refs. 25 and 27) and requires no 
modification for molecular ions, so only important steps and results are given in 
this appendix. 


Harmonic Oscillator Approximation 

There is general agreement (refs. 4, 5, 6, 7, 8, and 28) that the equilibrium con- 
figuration of Hg is an equilateral triangle with a proton at each vertex. Hence, 
belongs to point group Dg h . It follows by the methods of Wilson, Decius, and Cross 
(ref. 25) that the structure of the symmetry coordinate representation is 

T(s) = + E’ (Bl) 

where s stands for the three symmetry coordinates (ref. 25). The normal coordinates 
have the same structure of their representation as the symmetry coordinates do. There 
is one zero-order frequency associated with the normal coordinate of symmetry species 
Aj and another zero-order frequency associated with the pair of normal coordinates of 
symmetry species E'. In the following paragraphs, the zero-order frequencies of 
vibration will be found by setting up and solving the secular equation (ref. 25) in terms 
of the symmetry coordinates. 

The lengths of the three sides of the triangle with protons at the vertexes are called 
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inter nuclear distances. The changes in these internuclear distances from the equilibrium 
value r g are designated gp and g^, coresponding to sides 1, 2, and 3 of the 
triangle, respectively. With the aid of equation (Bl), the symmetry coordinates s are 
related to the changes in internuclear distances by the matrix equation 



S 1 


i/3 1 / 2 

1/3 1 / 2 

1/3 1 / 2 ' 


1 


S 21 

= 

2/6 1 / 2 

-1/6 1 / 2 

-1/6 1 / 2 


*2 


_ S 22_ 


0 

1/2 1 / 2 

-1/2 1 / 2 


kl 

where s* is of species AI and s 2 ^ 

and s ,, 

are of species E 

s were chosen to conform to Nielsen's notation in ref. 27. ) 



(B2) 


(The subscripts of 

The potential energy V and vibrational kinetic energy are also needed in 
setting up the secular equation. The symmetry of requires that, in the harmonic 
oscillator approximation, the potential energy (referred to the equilibrium value) be of 
the form 


2V’ - + <£3)+ ®2^1^2 + ^1^3 + ^2^3) 


(B3) 


where Bj and B 2 are constants. Combining equations (B2) and (B3) gives 


2V’ - B^ + s\ x + s| 2 )+ B 2 (s5 ' \ S 21 ‘ \ s 22) 


(B4) 


The kinetic energy is given by the matrix equation 

2E k = 1 Cf 1 I (B5) 

where G~* is the inverse of the G matrix (ref. 25) and the dots indicate derivatives 
with respect to time. Using the method in reference 25 results in 


2/m 

l/2m 

l/2m 

l/2m 

2/m 

l/2m 

l/2m 

l/2m 

2/m 


(B6) 


where m is the mass of a proton plus two-thirds the mass of an electron. From equa- 
tions (B2), (B5) and (B6), the following is obtained: 
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OT7 _ m ,1 -2 A 2 . 2 ^ 2 • 2 > 

2E k “ m (g S 1 + 3 S 21 + 3 s 22 


(B7) 


From equation (B4) and (B7), the secular equation is 


+ Bg — mx/3 
0 
0 


0 0 
B x - (B 2 /2) - 2mX/3 0 

0 B x - (B 2 /2) - 2mX/3 


= 0 


(B8) 


where the zero-order frequencies are given by X^ 2 / 27T. The solutions of equation (B8) 


are 


X l = 


3(Bj + B 2 ) 


m 


X 2 " 


3(B 1 - B 2 /2) 


2m 


(B9) 


the second solution occurring twice. 

The next problem is to obtain the L"* matrix (ref. 25), which is needed in the 
perturbation solution for the spectroscopic constants contained in a following section. 
Assume the molecular ion is in the x-y plane, with the axes passing through the center 
of mass as shown in the following sketch: 


y 
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where the protons are represented by circles and r is the equilibrium internuclear 

® th 

distance. Let Ax^, Ay^, and Az. be the components of the displacement of the j 
proton from equilibrium. Mass -weighted Cartesian displacement coordinates are defined 
by 


Vi = m 1/2 

tj 4 = m 1/2 Ax 2 , 
1/2 

Vrj= m AX 3 , 


77 2 = ml/2 Ayj, 
1/2 

% sm A ? 2 ’ 

1/2 

% 3 m ^ 


773 = m 1//2 Az j 
J7 6 Hm 1/2 Az 2 ► 

^9 s ml/2 Az 3 


(BIO) 


so that a matrix which appears later in the analysis will be orthogonal. For small dis- 
placements 


a 1 / 1 3!/ 2 1 \ 

= — ”1 + j ”2 + — ”4 - j ”5) 


m 


<£, = — 1 — , 

2 1/2 \ 2 


m 


t l/2 x 3 l/2 4 \ 

^l- 2 ^2 + — 


go = 


m 


1/2 


(Vq ~ V 5 ) 


(Bll) 


A set of coordinates such as s^, s^, s 22 in which the potential energy (eq. (B4)) 
and the vibrational kinetic energy (eq. (B7)) involve no cross terms (such as SjS 2 j) are 
always coordinates which differ from normal coordinates Q^, Q^, an< * Q22 on ^ a 
constant for each coordinate or do not differ at all. The constants can be found from the 
matrix equation 


(B12) 

by noting that L”* must be an orthogonal matrix. From equations (B2) and (Bll) and 
the proportionality between Sj and Q^, s 21 and Q 2 i> and s 2 2 and ^22’ e( l uation 
(B12) becomes 
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Qi 


Q 21 


CM 

CM 

<y 


q 4 


q 5 

11 

q 6 


s 


Q 8 


1 — 

£> 

CO 

1 



■S^/S 

0 

0 

3 1//2 /6 -1/2 

0 

s^/e 

1/2 

2^ 2 /6 

1/2 

0 

3 1//2 /3 0 

0 

-s^/e 

- 1/2 

-1/2 

CO 

CM^ 

> 

°? 

0 

0 3 1//2 /3 

0 

1/2 

-3 1/,2 /6 


01 

0 


**1 
^2 
% 

H 
% 

% 

Vrj 

^8 
(B13) 


where to Qg are normal coordinates for rotation and translation. The elements of 
L” ^ in equation (B13) that are pertinent to rotation and translation have been indicated by 
dashes since they are not relevant to finding spectroscopic constants. Thus the L - * 
matrix has been determined to the extent necessary in the following calculations. 


Evaluation of Potential Energy 


To obtain additional spectroscopic constants by the perturbation method of Nielsen 
(ref. 27), the cubic and quartic terms in the potential energy are needed in addition to 
the quadratic terms. All these terms are obtained by deducing the correct form of the 
potential energy and then selecting the constants to fit Conroy’s calculated potential 
(ref. 5). This is explained in more detail in the following paragraphs. 

The correct form for the potential energy is obtained by noting that the potential 
energy must be invariant for all symmetry operations permitted by point group 
even when Hg is vibrating. By application of group theory (refs. 13 and 25) this yields 


V = v e + ± 
e 2 


L l^l + ^2^21 + ^ 22 ) 


+ k lll Q l + k 122 Q l( Q 21 +Q L) +k llll Q l 

2 


2(2 2 \ (2 2 \ 

+ k n22 < ^l\ < ^21 + ^22/ + k 2222\^21 + ^22/ + hi & her order terms (B14) 


where V g is the potential energy at equilibrium, and k.^, k 122 > k im> k H22’ an< * 
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k 2222 are force constants. Expressing equation (B14) in terms of changes in bond 
lengths and omitting higher order terms gives 

V = V e + — - ( S\ + i\ + ^ 3 ) + + ^ 1^3 + ^2^2) + ( B 3 + B 4^(^l + ^2 + ^ 3 ) 


+ (6Bg - 3B 4 ) ^<^ 2^3 + ^ B 3(^l*^2 + ^1^2 + ^2*^3 + ^2^3 + ^3^1 + ^3^l) 


+ (Bg + Bg + B 7 )(^i + <^2 +, ^3) + ^5 + B 6 " 2B 7^(^ 1^2 + ^1^3 + ^2^1 

+ ^2^3 + ^3^1 + ^3^2) + ^ 6B 5 + 2B 7^('^1^2 + ^1^3 + ^2^3) 

+ (12 Bg - 3 B 6 )^ 2 ^ 2 -^ 3 +-^2^1^3 + ^ B15 ^ 


where Bg to B^ are given by 


B3 s m 


B 4 = m 


B 5 im 


B 6 ~ m 


B ? s m 


3/2 k lll 
27 

3/2 ^122 
27 

2 k llll 
81 

2 ^1122 
81 

2 16k 2222 
81 






(B16) 


Conroy’s computed Hg potential is shown in figure 2. He obtained an equilibrium 

internuclear distance r of 1. 68 Bohrs (8. 8900xl0” 9 cm) and a potential energy V at 

e 4 1 ^ 

equilibrium of -0. 357 Hartree (-7. 835x10 cm" ). Conroy very kindly furnished large 

graphs similar to figure 2 of this report and to figure 8 of reference 5. Potential 
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energies and internuclear distances were read from these two graphs by means of an 
optical comparitor. Points along lines J, Q, and S (fig. 2) with potentials of -0. 29 
Hartree (-6. 36xl0 4 cm -1 ) or less were used in fitting V (eq. (B15)) to Conroy’s poten- 
tial by minimizing the sum of the squares of the deviations. It was impossible to get a 
realistic fit by holding r and V at the values recommended by Conroy and varying 
Bj to B rj. Instead, the nine quantities r g , V e , and B^ to B ? were all varied simul- 
taneously. The resulting least-squares fit gave 


B ] 

B, 

B, 

t. 

B 4 

Bp 

B e 


B r 


1.66597 Bohrs (0. 88158xl0" 8 cm) 

-0.354771 Hartree (-7.7863xl0 4 cm" 1 ) 

0.2709 Hartree/Bohr 2 (2. 123xl0 21 cm'Vcm 2 ) 

-0. 11736 Hartree/Bohr 2 (-9. 198xl0 20 cm' 1 / cm 2 ) 
-0. 008914 Hartree/Bohr 3 (-1. 3203xl0 28 cm' 1 / cm 3 ) 
-0. 06664 Hartree/Bohr 3 (-9. 870xl0 28 cm" 1 / cm 3 ) 

0. 0016468 Hartree/Bohr 4 (4. 609xl0 35 cm'Vcm 4 ) 

0. 012395 Hartree/Bohr 4 (3. 469xl0 36 cm'Vcm 4 ) 

-0. 016949 Hartree/Bohr 4 (-4. 744xl0 38 cm'Vcm 4 ) 


(B17) 


Attempts to include points with larger potential energies in the fit were unsuccessful 
for some unknown reason. The reason was probably either lack of higher order terms 
in equation (B15) or inaccuracies in figure 2 of this report and figure 8 of reference 5. 

As might be expected, equation (B15) with constants from equation (B17) fit Conroy's 
potential (fig. 2) near the bottom of the potential well, as shown in figure 4, but not at 
higher potential energies. However, this was considered satisfactory because the fit is 
only used to calculate partition functions at temperatures low enough that most of the H^ 
ions occupy states near the bottom of the potential energy well. 


Perturbation Calculation and Final Results 

Nielsen (ref. 27) has carried out a second-order perturbation calculation to get the 
vibration -rotation energies of any polyatomic molecule except those with tetrahedral 
symmetry or internal rotation. In this section Nielsen's equations are applied to Hg to 
obtain its vibration-rotation energies and spectroscopic constants. Before doing this, 
it was necessary to correct a number of typographical errors in reference 27 in addition 
to the errors in reference 27 subsequently corrected by Nielsen (ref. 29). The necessary 
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force constants and equilibrium internuclear distance for Hg were obtained from the 
least-squares fit in the preceding section. 

It is customary in spectroscopic work to use term values rather than energy levels. 
A term value is defined as E/hc where E is the energy level, h is Planck’s constant, 
and c is the speed of light. Based on Nielsen’s equations and on the results of the sec- 
tion Harmonic Oscillator Approximation, the term values of Hg are 


E_ 

he 


he 


+ w lK + -J + w 2 (v 2 + 1 ) +x 



12 + X 21 J 


)(v 1 + — )(Vo + 1) 


+ x 22 (v 2 + 


i + j w * « b v - « 2 < b , 


°v> 


J 2 (J 


1) 2 D, 


- J(J + 1)k 2 D Jk - K 4 D k t 2C e Z 2 l 2 K 


(B18) 


where the quantum numbers v^, v 2 , 1 2 , k, and J take on the values 


v x = 0, 1, 2, 3, . . . 
v 2 = 0, 1, 2, 3, . . . 

l 2 =v 2’ v 2 ” 2 ’ v 2 " 4> ‘ ’ 
k = 0, 1, 2, 3, . . . 

J = K, K + 1, K + 2, . . . 




1, or 0 > 


(B19) 


and the other quantities are constants. The E Q /hc term in equation (B18) was not eval- 
uated by Nielsen and is probably very small. Consequently, it is neglected. The last 
term in equation (B18) accounts for the Coriolis splitting of the degenerate vibrational 
levels. Since the effect of Coriolis splitting on the partition function should be very 
small (ref. 30), the last term is neglected. 

Expressions for the spectroscopic constants in equation (B18) with the exception of 
£ 2 must all be found. From the section Harmonic Oscillator Approximation, 


U) l = 



2nc 


1/2 

2 

2 27TC 


(B20) 

(B21) 
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By using equation (B13) and Nielsen’s relations (ref. 27), the remaining spectroscopic 
constants are 


k ll 


32t7 4 w 2 c 3 


3k 


1111 


15k m 

877 2 c 2 w 2 


(B22) 


X 12 +X 21 


„„ 4 3 

3277 COjWgC 


2k 


3k lll k 122 


“122 


1122 


„ 2 2 2 2 2/, 2 2 

2n C 00^ 7T c (4ct>2 - w 


l) 


(B23) 


x 22 ■ 


90 4 2 3 
3277 WgC 


3k 


L 122 


2222 


, 2 2 2 2 2 \ 
477 c 0)^ 2(4^2 - 


(B24) 




k 2222 + 


k 122 


877 2 c 2 


( 4w 2 - w l) 


, 2 2 2 
477 C 0>2 

j( e ) 

zz . 


(B25) 


B v =B e -« 0 -(v 1 + ^« 1 -(v 2 + l) Q ! 2 

C v = C e-> / 0-( v l + ^l- (v 2 + 1) ’ / 2 


Dj =B 3 JUA 
d e 2 2 
V°1 w 2/ 


D t = -B 3 [ — + — 
Jk e l 2 2 

Vl w 2/ 


(B26) 

(B27) 

(B28) 

(B29) 


D =B 3 /-J- + 2 

x e l 2 2 

2co l w 2/ 


(B30) 


where, according to Nielsen (ref. 29), a Q and y Q are small and may be neglected and 
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(B31) 


B = 

-*>£?• 


a i = 


3B h 
e 


8„ 2 cl£>»i 


1 + 


hz 


1/2 


k lll 


o 2 2 2 
2ir c o>i 


a n = - 


B e h 


16» 2 cl £>» 2 


6 + 


1/2 


l 122 


2 2 2 
IT C Wj 




8 * 2 4z c 


V 1 


3B h 
e 




1 + 


:(e)l 1/2 k 

LZZ J 


ill! 


„ 2 2 2 
277 C Ci>i 


^2 


B e hk 122 


i« 4 3 

167T C 


r(e) 
zz . 


1/2 


w l w 2 


(B32) 


(B33) 


(B34) 


(B35) 


(B36) 


For calculating the partition function, the rotational -stretching constant p is 
required in Woolley’s approximation (ref. 12) rather than the constants Dj, Dj^, and 
D . These four quantities are related by the following equation (ref. 31): 


P 


^k_ 

he 


3(D J + D Jk + D k ) , 


4C 


2 

e 



2D J + D J K 
2C e B e 


(B37) 


Equations (B28) to (B30), (B34), and (B37) result in 



(B38) 
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The spectroscopic constants were evaluated from equations (B9), (B16), (B17), (B20) 
to (B25), (B31) to (B36), and (B38) and are 


= 3476 cm - * 

C = 21. 52 cm - * 
e 

cc >2 = 3601 cm - * 

q >2 = 1. 1848 cm - * 

= -70. 40 cm - * 

o »2 = 0. 1303 cm - * 

+ ~ -180. 43 cm - * 

yj = 0. 5924 cm - * 

Xg 2 = -76. 11 cm - * 

7^2 = 0- 8372 cm - * 

x, , = 23. 83 cm - * 

l T2 , 

p = 2.780X10 -5 °: 

B q = 43. 05 cm' 
e 



(B39) 


Values for o>j and Wg obtained by other investigators from calculated potential 
energies are given in table I. The agreement is not good, but the values (eq. (B39)) 
used in this report are probably the most reliable. It is probable that not one of the 
constants in equation (B39) is accurate to the number of significant digits given, because 
the Hg potential energy (ref. 5) was not evaluated for a sufficient number of isosceles or 
equilateral triangle configurations nor for any scalene triangle configurations. 

Since the ground electronic state of Hg is a *A^ state, it has an electronic degen- 
eracy g of 1. 

® . 

The energy D Q required to dissociate Hg in its rotationless ground vibrational 
state into two ground state H atoms and a proton (all infinitely separated) is obtained by 

subtracting the zero-point energy from the absolute value of V_ and is 72 707 centi- 

-1 e 
meters . 

From equations (B18) and (B39) the nondegenerate fundamental frequency oj \ is 

i 1 

3155 centimeters and the doubly degenerate fundamental frequency wjj is 3306 centi- 
meters”^ 1 . This completes the calculation of the spectroscopic constants and dissoci- 
ation energy of Hg. 
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APPENDIX C 

HIGH-TEMPERATURE APPROXIMATION FOR THE PARTITION FUNCTION OF Hj 

At temperatures of the order of 10 000° K and higher, the usual approximations 
(ref. 12) for the partition functions of polyatomic molecules and molecular ions do not, 
in general, give accurate results, nor does equation (B18) have sufficient terms to 
accurately predict all the energy levels of Hg that contribute significantly to the internal 
partition function. Instead, in this appendix the partition function of Hg is calculated 
classically, and corrections for the quantization of vibration are then applied. Due to 
necessary assumptions, the treatment is quite approximate. 


Analysis 

The determination of the internal partition function of Hg starts with equation (11-131) 
of Davidson (ref. 32). Davidson makes no distinction between bound states of a molecule 
and states in which the atoms have too much energy relative to one another to be bound. 
Hence, the partition function for all these three -atom states collectively is 




-5C/kT j3-r j3-r j3-r .3—- ,3 - 7 - j3rj- 

' d Pj d Pg d Pg d r^ d rg d rg 


(Cl) 


where q” t is the classical partition function including translation of the atoms, the r 
vectors all have an origin fixed in space, W is the volume of the container, and the 
number subscripts 1, 2, and 3 refer to the three atoms, respectively. The Hamiltonian 
3C is given by 


3C = 




v ft. r 2- r a ) 


(C2) 


If equation (C2) is substituted into equation (Cl), the momentum integrals can be 
evaluated with the result 


*ct = 


/27rmkT 


,9/2 


(C3) 
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where Z" is the configuration integral given by 


Z" 



-V/kT ,3— ,3- ,3- 

e a a Tg a rg 


(C4) 


Let rjg be the distance between atoms 1 and 2, r 13 the distance between atoms 
1 and 3, and r 23 the distance between atoms 2 and 3 (fig. 5). The crude assumption 
is made that V is the sum of three two-atom potential energies, irrespective of where 
the two electrons are; that is, 


V = U 12^ r 12^ + U 13^ r 13^ + U 23^ r 23^ 
Hill (ref. 26) has shown that exp(-U 12 AT) may be expressed as 

-u 12 /kT -uj 2 At -u| 2 At 


(C5) 


(C6) 


and similarly for U^ 3 and U 2 g. Here U?j is the effective potential energy between 
unbound atoms i and j and is given by 




W 


°ij< r ii> >0 


Oyfry) - kT in [1 - Ftry)] Uyfry) = 0 


(c 7) 


where the function F(r.j) has been defined by Hill (ref. 26). The quantity u£ is the 
effective potential energy between bound atoms i and j and is given by 


u i t i (r ii ) = 


U u( r «> > 0 

U y (r y ) -kTlnF^.) U^(r..) < 0 


i r 


ir i r 


(C8) 


Combining equations (C4) to (C6) results in 
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-< U f2 +U i3 +U l3> /kT -< u i2 +U 13* U 23>/ kT 

e + e 


-< U i2 +U !3 +U l3>/ kT -< u ! 2 +U !3 +U l3> /kT -< U 12 +U i3 +U l3> /kT 

+ G ~l“ G H - 6 

-< U ! 2 +U 13 +D 23> /kT -< u 12 +U 13 +U l3» /kT -< u 12 +U 13 +U 23 )/kT 

+ e + e + e 

d 3 ^ d 3 r 2 d 3 f 3 (C9) 


The first four terms in equation (C9) correspond to states in which at least one atom is 
not bound to any other atom and hence should be discarded. Integrals involving the next 
three terms are equal because of symmetry. Therefore, the appropriate configuration 
integral Z' for Hg is the sum of the last four terms in equation (C9) or 


Z' = 


Iff [-■ 


< u f2 +U i3 +U l3> /kT -< u i2 +u i3 +U i3> /kT ’ 


+ e 


'WW~W 


.3— .3— 

d r i d r 2 d r 3 


(CIO) 


12 ’ 


3_* 3— 3_ 

To evaluate equation (CIO), it is necessary to express d r^ d r 2 d r 3 in terms of 
The coordinate system for this transformation is shown in figure 5. 


r 13> and r 23- 


Proton 2 is the origin for polar coordinates r j 2 , 9^, cp^ of proton 1. An extension of 
r^ 2 serves as the z’ axis for polar coordinates r^ 3 , # 2 , <p 2 of proton 3 with proton 1 
as the origin. Obviously, 


d 3 Fj d 3 F 2 d 3 rg| = jd^rg r^ 2 dr^ 2 sin 6 ^ d^ d <p^ r^gdr 13 sin 0 2 d0 2 d(p 2 | (Cll) 


.3- 2 


The integrand in equation (CIO) is independent of r 2 , cp^, 6^, and <p 2 , so these six 
integrations can be carried out immediately. Integration over r 2 gives the volume W. 
Equation (Cll) becomes 


f J [ L |d 3 r 2 r2 2 dr 12 sine 1 d0 1 d<p 1 r 3 3 dr 13 sin0 2 de 2 d < p 2 | 

J(p 2 * / <p i ^ W 


= 8tt Wr 12 dr 12 r 13 dr 13 r 23 dr 23 


(C 12) 
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Combining equations (CIO), (Cll), and (C12) gives 


Z' = 8tt^W(3I 1 + I 2 ) 


(C13) 


where 


/ 


h s / r i2 e 


-U| 2 /kT /°°° -U{ 3 /kT /* r 12 +r 13 -u| 3 /kT 


L 


r 13 e 


L 


r 23 e "" dr 23 dr 13 dr 12 


12 -r 13 I 


(C14) 


and 


r -uf 2 /k T r 
h' / r 12* / 

Jo Jo 


-uf 0 /kT Z* 00 -ufg/kT / >r 12 +r 13 - U 23/ kT 


r 13 e 


jl2 +r li 
°^I r 12 _r 13 1 


r 23 e ” dr 23 dr 13 dr 12 


(Cl 5) 


The corresponding classical partition function including translation is 


q ct = 


27rmkT\ 9/ ^ 2 z , 
, h 2 ) 


(C16) 


(Compare with eq. (C3)). 

Equation (C16) was derived classically, so it requires two corrections: (1) it must 

be multiplied by the electronic degeneracy g g and (2) it must be divided by the 
symmetry number a. To obtain the internal partition function q’ , the resulting expres- 
sion must be divided by the translational partition function of the Hg ion; that is, 


_ 8l,2g e /27rmkT\ 3 
C ~3 3 / 2 a\ h 2 ) 


(3I : + I 2 ) 


(C17) 


The usual internal partition function q has the energy of the rotationless ground 
vibrational state as reference energy. Since q^. has infinitely separated H + H + H as 
the reference energy. 
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(C 18) 


-hcD /kT 

q c = q c e 

where D Q is multiplied by he because D Q is given in centimeters" 1 in appendix B. 

At low temperatures, equation (C18) is not accurate because the vibrations are 
quantized. As Davidson (ref. 32) has pointed out, this can be corrected for the 2 th 
normal vibration by multiplying the classical partition function by the ratio (R^ of the 
quantum -mechanical harmonic-oscillator partition function to the classical harmonic- 
oscillator partition function. 


-Uj/2 


®-l ~ 


n L e 


1 - e 


-u. 


where 


hca^ 

kT 


(C19) 


(C20) 


At high temperature, (R^ approaches 1. 

The final internal partition function q^+ is then given by 


q Hj " q c (R l (R 2 


87f2 g e 

3 3 / 2 <t 


/2mnkT\ 3 

V h 2 ) 


-hcD /kT „ 
(3I X + I 2 )e 0 (R X (R“ 


(C21) 


where (R 2 is squared because a> 2 belongs to a doubly degenerate vibration. 


Description of Digital Computer Program 

Equation (C21) was evaluated by means of a FORTRAN IV program for the Lewis 
IBM 7094 digital computer. 

The two-atom potential energy TJ ( r ) was obtained by taking one-third of Conroy's 
V(r^) for the equilateral triangular configuration (after adding 1. 00 Hartree 
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(2. 19475X10 5 cm -1 ) to Conroy's values, of course). These V’s are given on the line 

labeled Q in figure 2 and were read with an optical comparitor on large graphs similar 

to figure 2 of this report and to figure 8 of reference 5. These large graphs were very 

kindly furnished by Conroy. Thus U(r^) was obtained for internuclear distances from 

0. 4768x10“® to 1.9473x10”® centimeter. 

For larger and smaller internuclear distances, it was necessary to extrapolate 

U(r. •). For internuclear distances less than 0. 4768x10 ® centimeter, a Morse potential 

(ref. 33) passing through the points U(0. 4768x10”®) and U(0. 4937x10“®) and with the 

same potential well depth as U(r..) obtained from Conroy's data was used. For inter- 

nuclear distances greater than 1.9473x10” centimeter, another Morse potential passing 

through the points U(l. 8748x10”®) and 11(1.9473x10”®) and with the same potential 

well depth as U(ry) obtained from Conroy’s data was used. 

The triple integrals 1^ and 1 ^ were evaluated by numerical methods in the program. 

The values of q w + calculated by the program were mathematically accurate to within 
±lo 

0. 03 percent. 
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APPENDIX D 

COMPUTER PROGRAM FOR WKB ROTATIONAL ENERGIES OF DIATOMIC 
MOLECULES AND MOLECULAR IONS 

No complete table of vibration-rotation energies or term values of any diatomic 
molecule or molecular ion is available directly from an analysis of observed molecular 
bands. Such a complete table is necessary to calculate an accurate internal partition 
function for a diatomic molecule or molecular ion at high temperatures. Vibration- 
rotation energies may be calculated from theory if the potential energy U(r) of the 
molecule is known. For Hg the potential energy and vibrational energies are known 
from quantum theory (refs. 15 and 23). For many diatomic molecules, the most reliable 
potential energy curve that is obtainable is one calculated by the Rydberg-Klein-Rees 
(RKR) method (refs. 34, 35, and 36) which involves application of the WKB approximation 
(refs. 37 and 38) to observed molecular bands, the calculation of vibrational energies 
being an intermediate step. 

Since the vibration-rotation energy may be regarded as the sum of vibrational and 
rotational energies, a complete table of rotational energies (or term values) is generally 
the only missing data for calculating accurate internal partition functions. The program 
in this appendix calculates either a complete table of vibration-rotation term values or a 
complete table of rotational term values if a potential energy table is given as input. 

In the past, vibration -rotation energies of diatomic molecules have been calculated 
by at least three methods. Davidson (ref. 39) and Cooley (ref. 40) found a few energy 
eigenvalues of the appropriate second-order differential equation by supplying initial 
estimates of the vibration -rotation energy eigenvalues and iterating. Woolley, Scott, 
and Brickwedde (ref. 20) used the WKB connection formula to find noninteger values of 
the vibrational quantum number for given vibration -rotation term values and total angular 
momentum. 

For Hg three very sophisticated methods (refs. 15, 22, and 24) have been used to 
compute a few vibration -rotation energy eigenvalues to higher accuracy than for any 
diatomic molecule. 

For calculating partition functions, extremely accurate vibration -rotation energies 
are not required. Since in most cases the available potential energies are based on the 
RKR and WKB approximations, it is consistent to use the WKB connection formula to 
compute vibration -rotation energies. In this appendix the WKB connection formula is 
used to calculate vibration-rotation energy eigenvalues for any diatomic molecular ion or 
diatomic molecule to which the assumptions and restrictions listed in the next section 
apply. The described program has the unique advantages that it is unnecessary to know 
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beforehand if a state exists and that no initial estimates of vibration -rotation energies 
(or term values) are needed as input. 


Assumptions and Restrictions 

The assumptions and restrictions utilized in the program are listed herein for 
convenience; they are as follows: 

(1) Only *2 electronic states and electronic states obeying Hund's coupling cases (a) 
and (b) (ref. 33) are considered. (The ground state of Hg is a 2 2 state and obeys 
Hund's case (b)). 

(2) The extensions of the Born-Oppenheimer approximation (refs. 41 and 42) are 
assumed to be applicable. 

(3) The potential energy U(r) possesses only one well, as shown in figure 6, but it 
may have a hump (not shown). 

(4) The electronic state has no more than 100 vibrational energy eigenvalues (or 
term values) for its ground rotational state. 


Analysis 


In this section Hund’s coupling cases (a) and (b) are treated. When the total electron 
spin S and the component of orbital angular momentum along the internuclear axis A are 
zero, both coupling cases reduce to the case for ^2 states. The effective potentials 
for use with the WKE connection formula are given. 

Hund's coupling case (a). - In Hund's coupling case (a) the electronic motion 
including spin is coupled very strongly to the internuclear axis. From the relations 
given by Van Vleck (ref. 41) and Kronig (ref. 42) it can be shown that the differential 
equation for the radial wave function P(r) is approximately 



- W^(r) - B 


J(J + 1) - 2 A 2 


22ft + S(S + 1) 


+ 


47T<L • L) 



P(r) = 0 


pi) 


where ( ) indicates a time average over the electronic coordinates. 
The vibrational wave function ip i s defined by 
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if/ = r P(r) 


(D2) 


Combining equations (Dl) and (D2) gives 


3^ + 87T^m L _ w , (r) — 

dr 2 h 2 \ 87 r 2 /ur 2 L 


J(J + 1) - 2A Z - 22^ + S(S + 1) 


4tt 2 <L • L> 


i// = 0 (D3) 


The WKB connection formula (ref. 38) for equation (D3) is 


2(2 jihc) 


1/2 


r T ms 

•^min 


r max f 


/jr( v , J) - U^(r) - 




87r 2 jur 2 c 


J(J + 1) - 2A Z - 2SO 


+ S(S + 1) + 4 ’ 2 < z ■ l > 
h 2 . 


1/2 


dr = 


= H) h 


V = 0, 1, 2, 


(D4) 


where P is the vibration -rotation term value 


E F 

^-(v, J) = 

he 


P5) 


I and 


, v W a (r) 
u: (r) s — — 

a he 


(D6) 


and v is the vibrational quantum number. Equation (D4) can also be written 


2(2juhc) 1//2 /* rmax |>(v, J) - U”(r, J)] 1//2 dr = (v + |jh 
r min 


v = 0, 1, 2, . 


(D7) 
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where U’’ is an effective potential given by 

d 


U”(r,J) E U a (r) + 


hJ(J + 1) 

„ 2 2 
Sir lir e 


(D8) 


and 


u a ( r ) s U^(r) + 


87r 2 /ir 2 c 


-2A* 5 - 2SO + S(S + 1) + 


4n 2 (L • L) 


P9) 


The potential U (r) is the quantity obtained by applying the RKR approximation to a 

a 

molecular state obeying Hund's case (a). 

Hund's coupling cas e (jo). - In Hund's coupling case (b) only the electronic orbital 
motion is coupled very strongly to the internuclear axis. From the relations given by 
Kronig (ref. 42) and Mathews and Walker (ref. 38) it can be shown by a derivation 
similar to that for Hund's case (a) that the WKB connection formula is approximately 


2(2juhc) 


1/2 


f T max [>(v, K) - U£(r, K)] 1/2 dr = ^v + ±^h v = 0, 1, 2, 
' rc min 


PH)) 


where 


U£(r,K) =U b (r) + 


hK(K + 1) 

Sir 2 fir 2 c 


(DU) 


U b (r) s U b (r) 


87r 2 jur 2 c 


f-2A 2 + 4 * 2 < £ - L > 


(D12) 


The potential U b (r) is the quantity obtained by applying the RKR approximation to a 
molecular state obeying Hund's case (b). 

electronic state s. - The WKB connection formula for states may be obtained 
from equations (D7) to (D9) or (DIO) to (D12) and is as follows: 
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1/2 


r T max jy (v> J} _ u„ (r> J} j 1/2 dr = + 

r min 


v = 0, 1, 2, 


P 13) 


where 




U"(r, J) = U (r) + hJ(J + 


87r 2 Mr 2 c 


P14) 


and 


u » s V'(r) + <L ‘ 


2(J. r 2 ch 


P15) 


The potential U^Cr) is the quantity obtained by applying the RKR approximation to a 
1 S state. 

Generalized variables and rotation quantum number . - Inspection of equations p7), 
p8), plO), pll), pl3), and pl4) shows that the WKB connection formulas and the 
expressions for the effective potentials all have the same form, so that one computer 
program can be used for Hund's case (a) or (b) or for 1 E states provided the following 
correspondences (reading across) are kept in mind: 


Hund T s 

Hund’s 

1 s 

Generalized 

FORTRAN 

coupling 

coupling 

electronic 

variable or 

variable 

case (a) 

case (b) 

states 

quantum number 


U a 


U a 

U" 

UPP 

U a 

u b 

U a 

U 

U 

J 

K 

J 

j 

AJ 


In the remainder of this appendix, the generalized variables U" and U and the gener- 
alized rotation quantum number j are used instead of the corresponding quantities for 
particular coupling cases. 

Typical U and U" are shown in figure 6. In figure 6 and in the program, the 
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reference for energy is arbitrarily shifted to be U(r e ). If U” does not possess a well, 
there are no bound states. Even if U M does possess a small well, there will not be any 
bound state unless the left-hand side of equation (D7), (DIO), or (D13) (as appropriate) 
equals or exceeds h/2 when &= U”(r 04 ) (see fig. 6). 

The rotational term values f are found from 

f(v,j) + ^tv,j) - f(v, 0) (D16) 


Program Procedure 

The program used to calculate / or f was programmed in FORTRAN IV, IBM 
Version 13, for an IBM 7094 digital computer with the Lewis monitor system. Basically, 
the program interpolates a small table of U(r) given as input and obtains a large, 
equispaced table of U(r). For each value of j provided as input, the program computes 
a large, equispaced table of U"(r), determines if U"(r) possesses a well, and, if U" 
does, finds all Pty, j) by iteration of the generalized version of equations (D7), (D10), 
and (D13). These procedures are explained in detail in this section. 

Interpolation of U(r). - A set of points U(r) provided as input are represented by 
circles in figure 6. (In an actual problem many more points would normally be provided 
as input). The same interpolation formula is not used for all values of r. Instead, r is 
divided into regions 1, 2, 3, and 4 as shown in figure 6. In regions 1 and 4, the four- 
point interpolation formula of Lagrange is used. In region 2, a third-degree polynominal 
in r with zero slope at r is chosen to pass through U(r ) and the two points to the 
left of r g and is then used for interpolation. In region 3, a third-degree polynominal 
in r with zero slope at r is chosen to pass through U(r ) and the two points to the 
right of r g and is then used for interpolation. 

Interpolation of the points U(r) provided as input produces a large table of U(r) 
equispaced in r. The number of entries NMESH in this table is an input quantity. The 
first entry is the input point with smallest r, and the last entry is the input point with 
largest r. 

Procedure for each j_ value. - For each value of j provided as input, all solutions 
p (if there are any) are obtained and or f is printed out before going on to the next 
value of j. 

Calculation and examination of U"(r): The first step in this solution is to compute 
a value of U"(r) for each U(r) in the equispaced table by means of the generalized 
version of equations (D8), (Dll), and (D14). All curves U”(r) (fig. 6) possess two 
inflection points r Q1 and r^- By iterating on the second derivative of U”(r), the 
inflection point r Q ^ (fig. 6) is found. Examination of the first derivative of U”(r 0 j) 


38 


then reveals whether or not U”(r) has a well. If the first derivative is positive, there 
is a well; otherwise there is not. 

If there is a well, iterations on the first derivative of U"(r) provide r^g and r^^ 

(fig. 6). Setting equal to U"(r 0 ^) (see fig. 6) in the generalized version of 

equations (D7), (DIO), and (D13) gives the maximum value of v, which is generally a 

noninteger and is designated v . If v is negative, there are no eigenvalues 

max max 

If v max is 0, there is only one f (i.e., /=U"(r 04 )). If v max is positive, there are 

truncate (v ) + 1 eigenvalues P, where truncate ( ) means to drop any fraction in 
max 

*■ the argument. 

Table construction and search, and iteration: For v max that is positive, a table 
v search and one or two successive iteration procedures are used to find each eigenvalue 

f because the use of one iteration procedure alone is frequently inefficient or diverges 
in some case. Construction of the table, searching of the table, and iterating are 
explained in the following paragraphs. 

For use in the table searches, a table of v(£") is constructed from the generalized 
version of equations (D7), (DIO), and (D13) for equispaced & possessing values from 

U"(r n o) to U'^r,,.). The quantities r . (f) and r „„(/) are also included in the 

(Jo U4' min w max v 

table for use in the subsequent iterations. Although this table is adequate for table 
searches, v deviates substantially from being a linear function of so v and / are 
not appropriate variables for iteration. To overcome this, a Morse potential U^r) 

(ref. 33) is fit to each U”(r) so that the dissociation asymptote of the Morse potential 
has the value U"(rQ^). The vibrational term values ^.(v) of the Morse potential are 
used to obtain This is, in general, almost a linear relation, so and P are 

appropriate variables for iteration, ^ being given. The table for use in table searches 
therefore includes the five variables v, r min , r max > and It should be empha- 
sized that the use of is merely a mathematical device to make iterations converge 
more rapidly and in no way affects the accuracy of the resulting eigenvalues 

To find each eigenvalue 9" corresponding to an integer v, the table is searched to 
find the interval in the table wherein v falls. The value of ^ corresponding to the 
desired integer v is calculated and linear interpolation is used M times (M is an input 
) FORTRAN quantity) to get successively better values of P. After M linear interpola- 

tions, iteration is performed by using linear interpolation or linear extrapolation based 
, on the last two points calculated. If the new value of P so obtained falls outside the 

limits established in the table search, the new value of & is discarded, and five more 
linear interpolations are performed before attempting iteration by linear interpolation 
or linear extrapolation based on the last two points. Either iteration procedure is halted 
when the absolute value of the estimated error in !F is less than the input FORTRAN 
quantity TOL (the allowable error). 

Output options: After all eigenvalues ^(v, j) have been found for a given value of j, 
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the input FORTRAN quantities GV are subtracted from them . If all GV are zero, the 
output FORTRAN quantities FV will be the vibration-rotation term values P(v, j). If 
the GV are assigned the values f(?r, 0) determined from a previous run of the program, 
then according to equation (D16) the output quantities FV will be the rotational term 
values f(v,j). 

After writing (and, if specified, punching) the results for a given value of j, the 
processes are repeated for the next value of j . 

Organization of program and description of subprograms: A block diagram of the 
program and its subprograms is given in figure 7. SIMPS 1 is shown as a subprogram 
although at Lewis it is a built-in library function. Listings of the program and subpro- 
grams are given at the end of this appendix. 

The calls of the subprograms and the operations performed by them are as follows: 

The function LBINT has the call BINT(X, Y, Z, N) and performs four-point inter- 
polation by the method of Lagrange. X and Y are the independent and dependent 
variables, respectively, in a table of N entries. Z is the argument (value of X 
desired). BINT is the interpolated value of Y. 

The function LDER has the call DER(R) and is the first derivative of UPP with 
respect to internuclear distance R. 

The subroutine LDIFU has the call DIFU(RLOW, RHIGH, ROOT, $ STATEMENT NO. ) 
and finds the zero of UPP-FAA (FAA is a vibration-rotation term value in COMMON) 
which has R between RLOW and RHIGH. The value of R corresponding to the zero is 
ROOT. 

The function LEN has the call EN(R). EN is the value of the squareroot inside the 
integral in the generalized version of equations (D7), (DIO), and pl3). 

The subroutine LROOT has the call FIND(UUU1, UUU2,G,DELR, ROOT, $ STATE- 
MENT NO. ) and finds the zero of some specified function G(R) for R between UUU1 
and UUU2 by interval halving and linear interpolation. DELR is the spacing of the R 
table and is needed because G(R) is always derived from a table and is not continuous. 
The value of R corresponding to the zero is ROOT. 

The function LSECD has the call SECD(R) and is the second derivative of UPP with 
respect to R. 

The function LUP has the call UP(R). UP is the value of UPP at R obtained by 
linear interpolation of the equispaced UPP table. 

The subroutine LVW has the call VW(F, RMINH, RMINL, RMAXH, RMAXL, V, RMIN, 
UMAX, $ STATEMENT NO. ) and finds the vibrational quantum number V and corre- 
sponding near turning point RMIN and far turning point RMAX if the vibration-rotation 
term value F is given. It uses the generalized version of equations (D7), (DIO), and 
(D13). RMINL and RMINH are the values of R between which the search for RMIN is 
made. RMAXL and RMAXH are the values of R between which the search for RMAX is made. 


40 


The function SIMPS 1 has the call SIMPS 1(RMIN, RMAX, EN, L) and finds the integral 
of the function EN(R) between the limits R = RMIN and R = RMAX by a modified 
Simpson's rule. If an accuracy of 1 part in 10 cannot be achieved, 1 is added to L to 
flag the result. SIMPS1 is the value of the integral. This function was developed by 
T. E. Fessler and W. F. Ford of Lewis Research Center. 


Accuracy 

When TOL is set at its usual value of 0.25 centimeter - *, the accuracy of the com- 
puted ^(v, j) or f(v, j) is believed to be limited by the extended Born-Oppenheimer 
approximation, the WKB approximation, or insufficient points U(r) provided as input 
rather than being limited by inaccuracies in the program itself. For example, rotational 

term values for Hit (table V) computed by the program agreed within 0. 9 centimeter - * 

“ 1 , 
with rotational term values calculated from Wind's results (ref. 15). For the X 2g state 

of Hg the rotational term values agreed to within 7 centimeters - * with rotational term 

values calculated from wave numbers of lines measured by Herzberg and Howe (ref. 43). 

In general, slightly better accuracy is obtained if the program is used to compute 

f(v,j) rather than ^(v, j). If it is then desired to find ^"(v, j) by adding f(v,j) and 

^"(v, 0), values of 0) should be obtained from the same source as the potential U(r). 


Input 

Table VI shows the input variables required as they are to be punched on the data 

cards. The simplified sample problem given is Hg with only 30 potential energy points 

specified in addition to RE. The input FORTRAN variables are as follows: 

NL Number of members in UL and RL arrays; maximum of 100. 

UL Potential energy U for internuclear distances less than RE, measured 

from bottom of potential well (first value must be greater than dissoci- 
ation energy), cm - *. 

RL Internuclear distance r in monotonically increasing order; the RL array has 

a one-to-one correspondence with the UL array, A (1 A = 10"® cm). 

NR Number of members in UR and RR arrays; maximum of 440. 

UR Potential energy U for internuclear distances greater than RE, measured 

from bottom of potential well, cm - *. 
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RR Internuclear distance r in monotonically increasing order; the RR array 

has a one-to-one correspondence with the UR array (a suggested value 
for the last RR is 25), A (1 A = 10" 8 cm). 

NJ Number of members of AJ and NSTOP arrays; maximum of 45. 

AJ Generalized rotational quantum number j. 

NSTOP Number of mesh point where searches for zeros of first and second deriva- 

tives of effective potential with respect to internuclear distance are to 
stop; must be less than or equal to NMESH. The NSTOP array has a one- 
to-one correspondence with the AJ array. The purpose of this array of 
constants is to prevent searches at very large internuclear distances 
where the effective potential may be nearly constant, and consequently 
the four -point interpolation may cause erroneous values of the first and 
second derivatives of effective potential with respect to internuclear dis- 
tance. If there is no reason to suspect such errors or if there is no 
basis for estimating NSTOP, it should be set equal to NMESH. A value of 
NSTOP corresponding to AJ = 0 must always be set equal to NMESH if the 
potential possesses no hump. 

NV Number of members of GV array; maximum of 100. 

GV Vibrational term values measured from bottom of potential well, in mono- 

tonically increasing order. If the number of vibrational term values is 
unknown or if it is desired to have FV of output be vibration-rotation term 
values, put in 100 zeros and set NV = 100. If it is desired to have FV of 
output be rotational term values, put in values of GV equal to FV from a 
previous AJ = 0 run which had all GV values equal to zero, cm'l 

o o - 8 

RE Equilibrium internuclear distance, r g , A (1 A = 10 cm). 

RM Reduced mass of the two atoms of the diatomic molecule or molecular ion, 

including masses of electrons, ju, g. 

TOL Allowable absolute value of estimated error in FV of output (a suggested 

value is 0.25 cm - *), cm -1 . 

NZ Number of mesh points to be included in each increment in coarse search 

for negative second derivative of effective potential with respect to inter- 
nuclear distance (a suggested value is the nearest integer to NMESH/ 50). 

If no negative second derivative is found in the coarse search, a fine 
search is always made. 

NMESH Number of potential energy and effective potential mesh points; maximum 

5000 (a suggested value is 5000). 
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N 

NRITE 


NPUNCH 

M 

NAME 


Maximum number of iterations to find FV of output (a suggested value is 10). 

A value of 1 will cause NMESH values of internuclear distance RRR, poten- 
tial energy U, and effective potential UPP to be printed out. Zero will 
cause this output to be omitted. 

A value of 1 will cause output to be punched on cards in addition to being 
printed. Zero will cause the output shown in table VII to be printed only. 

Number of interpolations in iteration on FV before switching to interpolation 
or extrapolation based on the last two values (a suggested value is 5). 

Diatomic molecule and electronic state; this will appear on the printed out- 
put and on any cards punched out. 


Output 

The output for the sample problem is given in table VII. Each section of the output 
has been numbered to correspond to the following descriptions: 

(1) The first output of the program is a listing of the input data. 

(2) Output for AJ = 1, 30, and 41, respectively. There were 18, 7, and 1 rotational 
term values (eigenvalues) FV, respectively. The variables and their units are given in 
the next paragraph. 

(3) Output for AJ = 43. The effective potential had a well but no eigenvalue FV. 

(4) Output for AJ = 44. The effective potential did not have a well, so there were no 
eigenvalues FV. The slope of the effective potential at the inflection point r q-^ (fig. 6) 
was negative. 

(5) Appearance of output that should be obtained if the slope of the effective potential 
at the inflection point (fig. 6) was zero (this has never occurred). 

In the output with descriptions (2) to (5), the 14 unlabeled numbers that extend in a row to 
the far right are the flags M and L(l) to L(13), respectively. The flags L(l) to L(13) are 
described in the following paragraph. 

The variables unique to output are as follows: 

VV Vibrational quantum number, v. 

FV Eigenvalue, cm"\ If all GV's are zero, FV is the vibrational -rotation term 

value If the GV's are the vibrational term values, FV is the rotational 
term value f. 

ERROR Estimated error in FV, cm’\ 
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The quantities W, FV, and ERROR occur in triads along a row, each triad corre- 
sponding to an eigenvalue. 


V(l) 

F (1) 

V(100) 

F(100) 

AJ 

R4 

RRR(NSTOP) 

L(l) 

L(2) 

L(3) 

L(4) 

L(5) 

L(6) 

L(7) 

L(8) 

L(9) 


Value of v 

max 

Effective potential U”(r 04 ) measured from bottom of potential well 
(see fig. 6), cm“*. 

Value of v for bottom of effective potential well (always has value 
-0. 5). 

Effective potential U”(r o3 ) measured from bottom of potential well 
(see fig. 6), cm'''. 

Generalized rotation quantum number j . 

Internuclear distance r Q4 (see fig. 6), A (1 A = 10 cm). 

Internuclear distance r at which searches for zeros of first and second 
derivatives of effective potential with respect to internuclear distance 
stop, A (1 A = 10" 8 cm). For a valid solution R4 < RRR(NSTOP) 
unless (1) NSTOP = NMESH, (2) AJ is so small that there is little 
difference between U and UPP, and (3) UL has no hump. 

= 1 if search for negative second derivative of effective potential with 
respect to internuclear distance is unsuccessful. Normally, = 0. 

= 1 if search for inflection point r Q1 (fig. 6) is unsuccessful. 

Normally, = 0. 

= 1 if slope of effective potential at inflection point (fig. 6) is 

negative. Normally, = 0. 

= 1 if search for Tq 3 (fig. 6) is unsuccessful. Normally, = 0. 

= 1 if slope of effective potential at RRR(NSTOP) is positive or 0. This 
is satisfactory if (1) NSTOP = NMESH, (2) AJ is so small that there 
is little difference between U and UPP, and (3) UL has no hump. 
Normally, = 0. 

= 1 if search for r^ 4 (fig. 6) is unsuccessful. Normally, = 0. 

= 1 if search for r . for /(v m „, j) (see fig. 6) is unsuccessful. 

Ill XII Hid A 

Normally, = 0. 

= 1 if there are more than 100 eigenvalues FV for a given AJ, thereby 
exceeding the capacity of the program. Normally, = 0. 

= 1 if search for r min for & * f(v max > j ) is unsuccessful. Normally, 

= 0. 
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is unsuccessful. Normally, = 0. 


L(10) = 1 if search for r max 

L(ll) Positive if subprogram SIMPS1 is unable to integrate the generalized 

version of equation (D13) for /(v m „ v ,j) to an accuracy of 1 part in 
10° or better. Normally, = 0. 

L(12) Positive if subprogram SIMPS1 is unable to integrate the generalized 

version of equation (D13) for j) to an accuracy of 1 

part in 10° or better. Normally, = 0. 

L(13) = 1 if iteration to find FV does not converge in N iterations. Normally, 

= 0 . 

When NPUNCH = 1, values of VV, AJ, FV, and NAME will be punched on cards, in that 

order. 


Running Characteristics 

The time for the computer to execute solutions for 348 rotational term values with a 
tolerance of 0. 25 centimeter - '*' was 7 minutes. A case has not been found that required 
more than six iterations to converge to this accuracy. 

If the program runs improperly for no apparent reason, set NRITE = 1 and check the 
UPP array for spurious bumps, which may result from inaccurate or too small UR or 
UL arrays. 


45 



Listing of Program 


$ I BFTC MAIN 

C PROGRAM FOR WKB ROTATIONAL ENERGIES OF DIATOMIC MOLECULES AND 
C MOLECULAR IONS 
C 

EXTERNAL SECD,DER,EN 

DIMENSION RRR { 5000) t L ( 13) ,UL(100) .RLI100) ,URI440 
lit RR( 440) »AJ(45)»GVII00)»F(100) tVIlOOJ »F VUOO ) , VV< 100 ) ,NST0P<45> 
2, KM INC 1C0),RMAX(100) ,FFT(100> .FFFTUOOI >ERROR(100I 
COMMON/BLOCK l/UPP »RRR »RATI 0/BL0CK2/RM/ BLOCKS/ DELR* NMESH/BLQCK4/ 
1FAA/BL0CK5/L 

DOUBLE PRECISION UC5000) t UPPC5000) »8INT 
EQUIVALENCE (U( lit UP PCI) ) 

4 2 READ! 5,10CJNL,<UL!I i,I=l,NL) 

WRITE! 6, 1C1)NL,(UL(I ) ,I=1,NL> 

READ! 5,102HRL<IA),IA=1,NL) 

WRITE(6,103)CRLCIA),IA=1,NL) 

READ! 5tl00INR»(UR{IB)tIB=ltNR) 

WRITE! 6 t 105)NR y !UR( IB ) »IB=1 » NR) 

RE ADC 5>102) (RR (IC ) t IC =1 « NR) 

WRITE! 6* 1071 CRRCIC ),IC=1*NR) 

READ! 5.10CIN J,(AJ! IDI,ID=1,NJ) 

WRITE! 6y 1041 NJtCAJC ID),ID=1,NJ) 

READ! 5 1 106 1! NSTOP i I A A ) ,IAA=1,NJ> 

WRITE! 6* 108) (NSTOP (IAA) ,IAA=l f NJ) 

READ! 5, 10 0)NV,IGV! IE ) , IE=1 ,NV) 

WRITE! 6» 1 09) NVt ( G Vi IE),IE=1,NV) 

READ! 5,110)RE,RM,TOL,N2 f NMESH t N,NRITE,NPUNCH,M,NAME 
WRITE! 6, 111)RE,RM,T0L,NZ,NME SH , N, NRI TE .NPUNCH , M , NAME 
C END OF INPUT 

C CONSTRUCT INTERNUCLEAR DISTANCE RRR ARRAY 
DEL R=( RR! NR)-RL( 1 ) ) /FLOAT! NME SH-ll 
RAT I0=1«/DELR 
RRR! 1 )=RL C 1) 

DO 1 IG=2*NMESH 
1 RRR! IG)=RRR( IG-ll+DELR 
C PRINT OUT RRR ARRAY 



IF(NRITEM4,44, 

43 

43 

WRI TE( 6, 121 ) 



KB= 1 


30 3 

IFCKB-NMESH+fl) 

3CO,3CO,301 

300 

KC= KB+7 



GO TO 302 


301 

KC= NMESH 


302 

WRITEC 6, 1241KB, 

CRRR ( I ) , I =KB ,KC) 


IFIKC-NMESH) 304,44,44 

304 

KB= KB +8 



46 



II 


I * 


! : 


-v 


GO TO 303 

44 U{ 1 ) = UL( 1 ) 

C INTERPOLATE POTENTIAL UL IN REGION I 

NLAS1- INT ( 1.4(RL(NL)— RRRil) ) /DELR) 

00 2 I H=2»NLAS1 

2 U( IH) = BINT(RL,UL,RRR(IH) ,NL) 

C INTERPOLATE UL IN REGION 2 

NLAS2* INT(1.4(RE-RRR( 1))/DELR) 

NFIR2=NLAS1+1 

DENOM=!RL INL— 1 )— RE )**2*IRL (NL)— RE)**3— IRL! NL) — RE ) **2*C RL(NL-l)— RE) 
1**3 

A=( UL ( NL— 1 1* (RL (NL )— RE ) **3— UL! NL) * ( RL ( NL— 1 )— RE )**3>/DEN0M 
B=( UL ( NL ) *!RL( NL—ll— RE )**2— UL( NL— 1 )*(RL(NL)— RE)**2 ) /DENOM 
DO 3 I I=N F IR 2, NLA S2 

3 U( II) = (RRR(II )-RE)**2*(A+B*(RRR(II )-RE) ) 

C INTERPOLATE POTENTIAL UR IN REGION 3 

NFIR3=NLAS2+1 

NLAS3= INT ( 1. +(RR( 1 )— RRR (1) ) /OE LRI 

DENOM = (RR( D— RE 1**2* (RR ( 2) — RE ) **3— IRRI2)— RE)**2*(RR(1 )-RE) **3 
A=( UR ( 1 )*(RR ( 2 1— RE 1**3— UR (21 * (RRC L I— REl **3 ) /OENOM 
B— ( UR( 2 1* (RR ( 1 1— RE 1**2— UR (1>*(RR(2)— RE) **2 ) /OENOM 
00 4 I J=N FIR 3, NLA S3 

4 U( I J J = (RRR ( I J )— RE )**2* (A + B*(RRR(IJ)-REI) 

C INTERPOLATE UR IN REGION 4 

NFIR4=NLAS3+1 

00 5 IK-NFIR4.NMESH 

5 U( IK )=BINT IRR t UR »RRR C IK) »NR) 

DEBUG NLASlt NFIR2tNLAS2f NFIR3«NLAS3 ,NFIR4 
C PRINT OUT POTENTIAL U ARRAY 
IF(NRITE)46* 46,45 

45 WRITE! 6» 122) 

K B= 1 

305 IF(KB— NMESH+8) 306*306,307 

306 KC= KB 47 
GO TO 308 

307 KC= NMESH 

308 WRITE! 6, 125IKB , ( UII I , I =KB ,KC) 

IFIKC— NMESH) 309,46,46 

309 KB= KB 48 
GO TO 305 

46 REWIND 3 

WRITE! 3>!U(I )»I=1, NMESH) 

BACKSPACE 3 

C ENTER GENERAL IZEO-ROTATION— QUANTUM— NUMBER AJ LOOP 
00 6 J 1=1, NJ 
NSTP=N STOP! J I) 

00 7 IL=1 , 13 

7 LC IL )=C 

C CONSTRUCT EFFECTIVE POTENTIAL UPP ARRAY 
ANUM«2.799E— 23*AJ(JI)*(AJ( JI )4l,)/RM 
REAO! 3KUC I) ,1 = 1, NMESH) 

BACKSPACE 3 
00 8 IM«l, NMESH 

8 UPP! IM )=U( IM )4ANUM/RRR(I M)**2 
DEBUG AJ ( J I ) ,ANUM 

C PRINT OUT UPP ARRAY 


47 


4 


ll- 



I F| NR ITE)48, 48, 

47 

47 

WRITE! 6, 123) 



KB= 1 


310 

IFIKB— NMESH+8) 

311,311,312 

311 

KC= KB +7 



GO TO 213 


312 

KC= NMESH 


313 

WRITE I 6, 125) KB, 

IUPPCI ) , I =KB ,KC) 


IFIKC-NMESH) 314,48,48 

314 

KB= KB+8 



GO TO 310 



C FIND INTERNUCLEAR DISTANCE RN WHERE SECOND DERIVATIVE OF UPP WITH 
C RESPECT TO INTERNUCLEAR DISTANCE IS NEGATIVE 

48 TEMS=l.E+38 

DO 9 IN=NFIR3,NSTP,NZ 
TEM=SECDCRRR ( IN ) I 
IFC TEH 149,50,50 

49 RN-RRR(IN) 

GO TO 13 

50 IF! TEM— TEMSI 51,9,9 

51 T£MS=TEM 
KTEM=IN 

9 CONTINUE 

KT1=KTEM— NZ+1 
K T2=K T EM+NZ- 1 
DO 10 I0=KT1,KT2 
TEM = SECDI RRR CIO)) 

IF! TEM 152,10,10 
10 CONTINUE 
LI 11=1 
GO TO 11 
5 2 RN=RRR ( 10 ) 

C FIND INFLECTION POINT R1 

13 CALL F1ND(RE,RN,SEC0,DELR,R1,$12> 

DEBUG Rl.RN 

C FIND SLOPE OF UPP AT INFLECTION POINT R1 
SLOPE=DER (R 1 ) 

DEBUG SLOPE 
VI 100)=-. 5 

C TEST IF UPP HAS A WELL BY EXAMINING SIGN OF SLOPE 
IF! SLOPE) 14, 15,61 

61 IF(AJIJI) )62,62,16 

62 R3=RE 

GO TO 63 
12 L ( 2)=1 
GO TO 11 

14 LI 3)=1 
R4=0. 

GO TO 11 

15 VI 1)=- .5 
R4=0. 

F( 100 ) =UP IR 1 ) 

FI 1 )=F1 100 ) 

GO TO 41 

C FIND INTERNUCLEAR DISTANCE R3 CORRESPONDING TO BOTTOM OF UPP WELL 

16 CALL F INDIRE ,R 1»DER *DELR,R3,$17) 
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63 F( 100 ) =UP IR3 ) 

C FIND INTERNUCLEAR DISTANCE R4 CORRESPONDI NG TO TOP OF HUMP OF UPP 
SLQP=DER4 RRR (NSTP) ) 

DEBUG R3,FI100), SLOP 
I FI SL0P)19, 18 t 18 

17 L I 4 ) = 1 
GO TO 11 

18 LI 5 ) = 1 
R4=RRR INSTP ) 

GO TO 53 

19 CALL F IND4R 1 f RRR INSTP ) ,DER ,DELR *R4 9 $20) 

53 FI 1 )— UP! R 4 ) 

DEBUG R4, F(l) 

FAA=F I 1) 

C FIND NEAR TURNING POINT RMINI1) FOR TERM VALUE EQUAL TO UPP AT R4 
CALL OIFUCRRRI 1) ,R3,RMINI 1) ,$21) 

RMAXI 1 )=R 4 

C FIND VIBRATION NUMBER V FOR TERM VALUE EQUAL TO UPP AT R4 

VI 1)=SQRTI3.620E+21*RM)*SIMPS14RMIN41) 9 RMAXI l ) »EN,LI11) >-.5 
DEBUG VI 1) 9 RMIN(1) 

I FI VI 1 >-99.999)40,409 22 

20 L ( 6 )= 1 
GO TO 11 

21 L( 71=1 
GO TO 11 

22 L ( 8 )- 1 
GO TO 11 

C TEST IF THERE ARE ANY EIGENVALUES I NONE IF VII) IS NEGATIVE) 

40 IF! VI 1 ) )4 1 9 23, 23 

41 WRITE I fc, 116) VI 1),F(1) , VI 100) ,F 1100) 

GO TO 11 

C CONSTRUCT TABLE OF TERM VALUES F, VIBRATION NUMBERS V, NEAR TURNING 
C POINTS RM IN, AND FAR TURNING POINTS RMAX 

23 NN= INT 1 VI 1)1 
IFINN- 10)25,26,26 

25 NN= 10 

26 DEL F=I FI 1 )-F 1 100) ) /FLOAT INN) 

IQ-NN + 1 

DO 27 IP = 2 , I Q 

27 F( IP ) = F( IP— 1 )-DELF 
DEBUGIF4IR),IR=1,IQ) 

VI IQ ) . 5 

RM INI IQ )=R 3 
RMAXI IQ )=R3 
DO 28 NM=2,NN 

CALL VVVIFINM) ,R3,RMINCNM— 1) ,RMAXINM-1) ,R3,VINM) ,RMINC NM) , 

1RMAXI NM),$11) 

DEBUG FINM) , VINM) ,RMININM) , RMAX I NM) 

28 CONTINUE 

C FIT MORSE POTENTIAL 
DI SE=F 1 1 )— FI 100) 

FAA=Fl 100)*DISE*.5 
CALL D IFU IR3,R4,R 5,$2 1 ) 

OM£G= 1*30 IE- 11*SQRT 1DI SE /RM)/IR5-R3) 

OM EGX=40MEG* I VI 1 ) + • 5)— DI SE ) /I VI 1 ) + . 5) **2 
DEBUG DIS£,FAA ,R5,0MEG ,OMEGX 
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C ADO NORSE TERN VALUE FFT TO TABLE 
00 700 NQ = 1, IQ 

700 FFT1NQ )=! VINQ )+. 5)* IONEG-ONEGX* ( V !NQ)+. 5) I 
DEBUG !FFTII),I=1,IQ) 

C CONSTRUCT TABLE OF DESIRED VALUES OF VIBRATION QUANTUN NUNBER VV AND 
C NORSE TERN VALUE FFFT 
32 VV! 1)=C. 

FFFT! 1 )=0MEG*.5— 0NEGX*.25 
IQ=INT (VC lll+I 
IF( 10-1)70,70,85 
85 DO 35 NS=2, IQ 

VV!NS)=VV!NS-1) + 1. 

35 FFFT!NS)=! WINS )+. 5)* ! ONEG-ONEG X*C V V I NS) *. 5 ) ) 

DEBUG! VV! I), 1=1, IQ) 

DEBUG! FFFT!I),I=1,IQ) 

C ENTER LOOP THAT SOLVES FOR ALL EIGENVALUES FV FOR GIVEN UPP 

70 DO 97 IZ=1,IQ 
NRR=0 

C SEARCH Ft V, FFT TABLE 
K=2 

71 IFtVIK)— VV!IZ))73,72,72 

72 K=K+1 

GO TO 71 

73 FFT 1=FFT! K ) 

FFT2=FFT! K— 1 ) 

F1=F( K ) 

F2=F( K— 1) 

FFTNEW=FFT1 
FNEW=F 1 

C ENTER ITERATION LOOP FOR EIGENVALUES FV 

701 FOL D=FNEW 

FN£W=F J+!F2-F1)*!FFFT! IZ)— FFT1 ) / CFFT2— FFT1 ) 

77 NRR=NRR+1 

IF! NRR— N ) 79, 79,78 

78 LI 13 ) = 1 
GO TO <7 

79 CALL VVV( FNEW,RNIN!K) , RNI N! K— 1 ) , RNAX 1 K— 1 ) , RNAX! K) , V3 , RNNIN, 
1RNNAX, $11) 

FFTOLD=FFTNEW 

FFTNEH=! V3+. 5 )* ! ONEG-ONEGX* ! V3+.5I ) 

SLOPP = !FOLD-FNEH) /1FFT0LD-FF TNE W) 

ERROR! IZ ) = ! FFTNEW— FFF T(IZ) )*SLOPP 

DEBUG FNEW.FFFTUZ), FFTNEW. S LOPP .ERROR 1 1 Z) ,FFT 1,FFT2,F1,F2,N 
IF! ABS! ERROR ! IZ ) )-TOL) 702,702,74 

702 FV! IZ )=FNEW-GVI IZ) 

GO TO S7 

74 IF! NRR— N )80,7C3,703 

80 IF! ERROR! IZ))81,81,82 

81 FFT1=FFTNEW 
F1=FN EW 

GO TO 701 

82 FFT2=FFTNEW 
F2=FN EW 

GO TO 701 

703 FOL D=FNEW 

FNEW=FOLD+SLOPP* ! FFF T ! IZ)— FFTNEW) 
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IFCFNEW-F(K ) ) 704, 705,705 
705 IF! FN£W-F(K-1))77,77,704 
704 M=M + 5 

FNEW-FOLD 
GO TO 80 
97 CONTINUE 

C PRINT OUT AND PUNCH ANSWERS 

WRITE ( 6,1X71 (VVCI)fFVfl) , ERR OR !I),I~1,IQ) 

IFINPUNCH >201,41,201 
201 00 200 NA = 1,IQ 

200 PUNCH 120 , VV CNA ) ,A Jl JI ) ,F V! NA) , NAME 
GO TO 41 

11 WRITE! 6, 118JAJCJIJ , R4,RRR( NSTPI ,N, I LII) ,1=1,13) 

6 CONTINUE 
GO TO 42 

100 FORMAT ( 1 10,7F10.0/( 8F 10* 0) ) 

101 FORMAT ( 1H 1, 3H NL,I3,3H UL, 1P9E13.6/ Cl X.9E13.6 I) 

102 FORMAT ! 8F10. 0) 

103 FORMAT 1 1H0,3H RL, 1P9E 13. 6/C 1 X, 9E13. 6) ) 

104 FORMAT C 1H0,3H NJ,I3,3H A J , 1 P9E13. 6/ (1 X .9E13.6 1 ) 

105 FORMAT ( 1H0,3H NR,I3,3H UR , 1P9E13. 6/ Cl X ,9E13.6 >) 

106 FORMAT (1615) 

107 FORMAT C 1H0,3H RR , 1P9E 13. 6/ C 1 X, 9E13. 6) I 

108 FORMAT ( 1H0, 6H NSTOP, 241 5/ IX, 211 5) 

109 FORMAT C 1H0 , 3H NV,I3,3H G V , 1P9E 13. 6/ Cl X ,9E13.6 )) 

110 FORMAT !3E13.6,6I5,A6) 

111 FORMAT C 1H0, 3H RE, 1PE13.6,3H RM,E13.6,4H T0L,E13.6,3H NZ,I5,6H NMES 
1H, I5,2H N,I5,6H NRITE,I5,7H NPUNCH,I5,2H M, 15 ,2X , A6// I 

116 FORMAT C 1H0, 5H VC1),F9.5,5H F (1 ) ,F1 0.2 ,9X,7H V C 100 ) ,F9. 5 ,7H FC100I* 
IF 10 .21 

117 FORMAT C 1H0, 3 C 3H VV,F4.0,3H FV,F10.2,6H ERROR, F9 .2, 5X1/ C IX , 3C 3H VV, 
1F4.0, 3H FV,F10.2,6H ERROR ,F 9.2 ,5X1)1 

118 FORMAT ( 1H0, 3H AJ,F1C.5,3H R4, 1 PE14. 7 ,1 1H RRR C NSTOPI ,E14.7 , 1415// I 

120 F0RMATIF4.0, 2X,F 10.5, 2X, FI 0.2,2 X,A6> 

121 FORMAT C 1H0,4H RRR) 

122 FORMAT ! 1H0,2H U) 

123 FORMAT! 1H0,4H UPP) 

124 FORMAT! 1H , I 4, 1P8E 14. 7) 

125 FORMAT! 1H , 14, 1P8014.7) 

END 
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I I I ill II 


iimuii 


$ IBFTC LBINT 

SUBPROGRAM FOR 4 POINT INTERPOLATION 

DOUBLE PRECISION FUNCTION BINT(X,Y,Z,N) 

DIMENSION X( 440) » YI440) 

DOUBLE PRECISION Q1»«2,Q3,Q4,Q12,Q34,P21.P32,P42 ,P31,P41,P43 
Ml=l 
M2=N 

5 I FI M2— Ml— 1 } 7 * 7 » 6 

6 MT£M=M1+INT( .5* (FLOAT I M2-M1 )+. 0001 ) I 
IFl 2— XIMT EM ) ) 8» 1 » 9 

8 M2=MTEM 
GO TO 5 

9 M1=MTEM 
GO TO 5 

1 BINT=Y(MTEMI 
RETURN 

7 IF(M1-1)2,2,3 

2 1=2 
GO TO 4 

3 IFCM2— N)ll,10,10 

10 I=N— 2 
GO TO 4 

11 I=M1 

4 Q1=Z-X(I-1) 

Q2=Z— XI I) 

Q3=Z— X ( 1+ 1 ) 

Q4=Z— X ( 1+2 ) 

Q 12=Q 1*Q2 
Q34=Q3*Q4 
P21=XI I)-X(I-1I 
P32=X( I*1)-X(I) 

P42=X( I+2I-XC IJ 
P31=X( I+n-XII-l) 

P41=XI 1+21— XI I— II 
P43=X( 1*2 I-XII+l) 

BINT=Q34*IQ1* Y 1 1 ) /( P32*P42)— Q2* Y(I — 1)/(P31*P41) )/ P21+Q12*! Q3*Y 1 1+2 
l)/( P41+P42)— Q4*Y( I+ll/I P31*P32) ) /P43 
RETURN 
END 
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$ I BFTC LDER 


C SUBPROGRAM FOR FIRST DERIVATIVE OF UPP 
C 

FUNCTION DER(R) 

COMMON /BLOCK 1/ UPP ,RRR,RATI 0/BL0CK3 /DELR, NMESH 
DIMENSION RRR( 5000) 

DOUBLE PRECISION UPP I 5000) 

FL0M = 1 • + ( R-RRR ID )*RATIO 
1= INT ( FLOM) 

IF( I— NMESH ) 1 9 2 » 2 

2 DER=( UPP ( NMESH )-UPP(NME SH— 1 J ) /DELR 
RETURN 

l DER=( UPP ( I + D— UPPII)) /DELR 
RETURN 
END 


$ I BF TC LDIFU 


C SUBPROGRAM TO FIND ROOT OF UPP-FAA 
C 

SUBROUTINE D IFU CRLQW , RHIGH .ROOT,*) 

COMMON /BL0CK4/FAA /BLOC Kl/UPP , RRR ,R ATI 0/BL0CK3 /DELR. NMESH 
DIMENSION RRR ( 5000) «NU ( 2 ) *F (2) 

DOUBLE PRECISION UPP (5000) 

NU( 1) = INT (1* +( RLQ W— RRR (I)) *RATI 0) 

NU( 2 ) = INT ( 1«+{RHIGH-RRR( 1) )*RATI0)+1 
I F( NU { 2)-NMESH)21,21,20 

20 NU( 2 ) =NME SH 

21 DO 100 1=1,2 
NV=NU( n 

100 F( I ) = UPP ( N V)— FAA 

3 IFiFI 11)4,22,5 

4 IF( F( 21)13,14,23 

13 RETURN 1 

14 N 1=NU ( 2) 

15 RQOT=RRR( N1 ) 

RETURN 

22 N1=NU( 1) 

GO TO 15 

5 IF(F( 2)118,14,13 
18 NT EM 1=NU( 1) 

T EM 2= F ( 1 ) 

NU( 1 ) =NU( 2 1 
F( 1 ) = F ( 2 ) 

NU( 2 ) =NTEM 1 
F( 2) = TEM2 

23 NV = INT(( FLOA T( NU( 1 ) +NU ( 2 ) ) + . 0001 ) /2. ) 

FUNC= UPP < N V ) — FAA 

9 i F( FUN 010,16,11 

10 NU( 1 ) =NV 
F( n = FUNC 
GO TO 24 

16 RQOT=RRR ( N V) 

RETURN 

11 NU( 2) =NV 
F( 2 ) = FUNC 

24 IF ( IABSiNUC 1 l-NUI 2) J-l ) 25 ,25 ,23 

25 N1=NU( 1) 

N2=NU (2) 

ROOT=( RRR (N1 )*F ( 2 J-RRR ( N2 ) *F ( 1 I ) / (F (2I-FC1) ) 

RETURN 

END 
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now o n 


$ I BFTC LEN 


SUBPROGRAM FOR SIMPS1 INTEGRAND 

FUNCTION ENI R) 

COMMON /BLOCK 4/F 
RR=R 

A*F~UPIRRI 
IF( A) 1 * It 2 

1 EN=0. 

RETURN 

2 EN = SQRT( A I 
RETURN 
END 


BFTC LROOT 

SUBPROGRAM TO FIND ROOT OF GCR) 

SUBROUTINE F INDt UUU1 , UUU2 ,G ,DE LR,RGOT,*l 
DIMENSION U4 2) t F ( 2) 

Ul 11*UUU1 
Ul 2)=UUU2 

JMAX=1 + INT( ALOGI I U € 2 ) — U ( 1 1 I /DELR) /*6931l 
1=0 

1 1 = 1+1 
V=U(I ) 

100 FUNC=Gf V ) 

I FI 2- I )9 f 2 f 2 

2 FI I ) = FUNC 
IF! 2-1 )3, 3 f I 

3 1=1+1 

IFI Ft 1))4, 15, 5 

4 IFI Ft 2 ) ) 13, 14, 7 

13 RETURN 1 

15 R00T*U(1J 
RETURN 

14 R00T=UC 2) 

RETURN 

5 IF(F(2))ie,14,13 
18 T EM 1=U ( 1 ) 

TEM2=F ( 1) 

U( 1 ) = U( 2) 

Ft 1 )=F( 2) 

Ul 2 ) = T EM 1 
FI2) = TEM2 

7 J=0 

8 J=J+1 

V=IUI 1 )+U ( 2) )/2. 

GO TO 100 

9 1 FI FUN 010,16,11 

10 Ul 1 ) = V 

F( 1 ) = FUNC 
GO TO 12 

11 U(2) = V 

FI 2 )=FUNC 

12 IFIJMAX-J 117,17,8 

16 ROOT=V 
RETURN 

17 ROOT* | Ul 1 )*F 12 )“UI 2)*F II ) 1/ IF 12)— Fll) I 
RETURN 

END 
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no & no 


SIBFTC LSECD 


SUBPROGRAM FOR SECOND DERIVATIVE OF UPP 
FUNCTION SECDCR) 

COMMON /BLOCK I /UPP , RRR , RAT I 0/BL0CK3 /DELR, NMESH 
DIMENSION RRR ( 5000 
DOUBLE PRECISION UPP (5000) 

FLOM= 1 •♦( R— RRR ( I ) )*RA TIO 
I=INT( FLOM) 

I F( 1-111, 1,2 

1 SECDM UPPC3I-2.*UPP(2) + UPP(1)I /DELR**2 
RETURN 

2 IF( I— NMESH 14,3, 3 

3 SECD= ( UPP (NMESH)— 2«* UPP ( NME SH— l )♦ UPP ( NMESH— 2) ) /DELR**2 
RETURN 

4 SECD=( UPP ( 14 1)-2.*UPP ( I )*UPP(I-1) ) /DELR**2 
RETURN 

END 


IBFTC LUP 

SUBPROGRAM TO LINEARLY INTERPOLATE UPP 
FUNCTION UP (R ) 

COMMON /BL0CK1/UPP, RRR, RATIO 
DIMENSION RRR (5000) 

DOUBLE PRECISION UPP(5000) 

FL0M=1 .+( R-RRR (1) )*RATIO 
I^INTC FLOM) 

UP^UPP (I)+(UPP( I+1)-UPP(II )*(FLOM-FLQAT(I ) ) 

RETURN 

END 


SIBFTC LVVV 

C SUBPROGRAM TO FIND VIBRATION NUMBER (OR QUANTUM NUMBER) FOR GIVEN TERM 
C VALUE F 
C 

SUBROUTINE VVV( F,RMINH,RMI NL, RMAXH ,RMAXL, V, RMIN, RMAX,*) 

COMMON /BL0CK4/FAA /BLOC K2/RM/BL0CK5/L 
DIMENSION L( 13) 

EXTERNAL EN 
FAA=F 

CALL D IFUCRM INL , RMI NH ,RMI N , $29) 

CALL DIFUIRMAXL, RMAXH, RMAX, $30) 

V=SQRT (3« 620E*21*RM)* SIMPS1 (RMI N,RMAX,EN,L«12 ) )-,5 
RETURN 
2S L ( 9 )= 1 
RETURN 1 
30 L ( 10)=1 
RETURN 1 
END 
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TABLE I. - COMPARISON OF CONSTANTS FOR THE EQUILATERAL TRIANGLE 


CONFIGURATION OF Hg OBTAINED IN VARIOUS CALCULATIONS 


Investigators 

Equilibrium 

internuclear 

distance, 

r e’ 

cm 

Equilibrium 

potential 

energy, 

V hc ’ 

cm - * 

Zero-order 

vibrational 

constant , 

-1 

cm 

Refer- 

ence 

w i 

w 2 

Hirschfelder 

9. 47x10" 9 

-6. 439X10 4 

1550 

1100 

3 

Ellison, Huff, and Patel 

9.31 

-7. 832 

3454 

2326 

8, 44 

Christoffersen 

8.771 

-7. 3006 

3354 

2790 

4 

Conroy 

8.89 

-7.835 

(a) 

(a) 

5 

Hoyland 

8.89 

-5.983 

(a) 

(a) 

6 

Pearson, Poshusta, and Browne 

8.78 

-6.990 

3610 

4440 

7 

This report 

8.816 

-7.786 

3476 

3601 

— 


a Not given. 

k Using least-squares fit to Conroy's potential energy (ref. 5). 



TABLE H. - PARTITION FUNCTIONS OF Hg AND Hg 


Temperature, 

Partition 

Partition 

Temperature, 

Partition 

Partition 

T, 

function of 

function of 

T, 

function of 

function of 

°K 

H 3> 

H+ 

°K 

H+ 

H 2> 





+ CO 

cjp 


298.15 

4. 8065 

7.5178 

6 500 

5. 9064X10 3 

6. 3778xl0 2 

300 

4. 8506 

7. 5627 

7 000 

8.2440 

7.4813 

400 

7.4266 

1. 0001X10 1 

7 500 

1. 1374X10 4 

8,6603 

500 

1. 0359X10 1 

1.2472 

8 000 

1.5520 

9.9048 

600 

1.3623 

1.4998 

8 500 

2.0951 

1. 1205X10 3 

700 

1.7220X10 1 

1. 7606x10* 

9 000 

2.7988X10 4 

1.2552X10 3 

800 

2.1167 

2.0319 

9 500 

3.7011 

1.3937 

900 

2.5503 

2.3159 

10 000 

4. 8460 

1. 5351 

1000 

3.0275 

2.6143 

11 000 

8.0735 

1.8238 

1100 

3. 5544 

2.9285 

12 000 

1.2969X10 5 

2.1165 

1200 

4. 1372X10 1 

3.2597X10 1 

13 000 

2. 0129X10 5 

2. 4092X10 3 

1300 

4.7830 

3.6088 

14 000 

3. 0243 

2.6991 

1400 

5.4991 

3.9768 

15 000 

4. 4063 

2.9839 

1500 

6.2936 

4. 3643 

16 000 

6.2364 

3.2623 

1600 

7. 1746 

4.7721 

17 000 

8.5881 

3. 5332 

1700 

8. 1509X10 1 

5.2010X10 1 

18 000 

1. 1525X10 6 

3.7959X10 3 

1800 

9.2319 

5.6516 

19 000 

1.5091 

4. 0501 

1900 

1.0427X10 2 

6. 1245 

20 000 

1.9310 

4.2955 

2000 

1. 1748 

6.6205 

22 000 

2.9635 

4.7604 

2200 

1.4809 

7.6846 

24 000 

4.2018 

5. 1916 

2400 

1. 8511X10 2 

8.8495X10 1 

26 000 

5. 5432X10 6 

5. 5908X10 3 

2600 

2.2965 

1.0121X10 2 

28 000 

6.8387 

5.9605 

2800 

2.8294 

1. 1507 

30 000 

8.2173 

6.3031 

3000 

3.4637 

1.3014 

32 000 

9.6610 

6. 6208 

3200 

4.2149 

1. 4648 

34 000 

1. 1154X10 7 

6.9160 

3400 

5. 1005X10 2 

1.6417X10 2 

36 000 

1.2682X10 7 

7. 1906X10 3 

3600 

6.1401 

1.8327 

38 000 

1.4233 

7. 4465 

3800 

7.3555 

2.0385 

40 000 

1.5797 

7.6855 

4000 

8.7712 

2.2595 

43 000 

1.8147 

8.0154 

4300 

1. 1331X10 3 

2.6207 

46 000 

2.0485 

8. 3149 

4600 

1. 4509X10 3 

3.0184X10 2 

49 000 

2.2790X10 7 

8. 5879X10 3 

4900 

1.8428 

3.4530 

52 000 

2.5052 

8. 8376 

5200 

2.3235 

3.9247 

56 000 

2.7985 

9. 1390 

5600 

3.1320 

4.6104 




6000 

4. 1768 

5.3590 
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TABLE m. - THERMODYNAMIC PROPERTIES OF Hg 


T, 

°K 

R 

h t -h 0 

RT 

s° 

IT 

g t- h o 

RT 

fr 

RT 

_ G T 

RT 

296.15 

4.0176 

3.9745 

17.7831 

13.8086 

441.1777 

— 423.3946 

300 

4.0177 

3.9748 

17.8080 

13.8332 

# 38.4819 

- 420.6739 

400 

4.0264 

3.9865 

18.9648 

14.9783 

329.8668 

- 310.9020 

500 

4.0517 

3.9965 

19.8655 

15.8689 

264.7008 

- 244.8353 

600 

4.1113 

4.0101 

20.6088 

16.5987 

221.2637 

- 200.6548 

700 

4.2132 

4.0314 

21.2497 

17.2184 

190.2487 

- 168.9990 

800 

4.3535 

4.0625 

21.8211 

17.7586 

167.00 27 

- 145.1816 

900 

4.5216 

4.1040 

22.3433 

18.2394 

148.9397 

- 126.5963 

1000 

4.7057 

4.1548 

22.8292 

18.6743 

134.5070 

- 111.6778 

1100 

4.8959 

4.2136 

23.2866 

19.0730 

122.7155 

- 99.4289 

1200 

5.0849 

4.2783 

23.7207 

19.4424 

112.9051 

- 89.1844 

1300 

5.2680 

4.3474 

24.1350 

19.7875 

104.6183 

- 80.4833 

1400 

5.4424 

4.4195 

24.5318 

20.1123 

97.5281 

- 72.9963 

1500 

5.6067 

4.4932 

24.9130 

20.4198 

91.3946 

- 66.4817 

1600 

5.7606 

4.5677 

25.2798 

20.7121 

86.0378 

- 60.7580 

1700 

5.9042 

4.6421 

25.63 34 

20.9913 

81.3198 

- 55.6864 

1800 

6.0380 

4.7160 

25.9747 

21.2587 

77.1338 

- 51.1591 

1900 

6. 1626 

4.7889 

26.3045 

21.5156 

73.3953 

- 47.0907 

2000 

6.2789 

4.8605 

26.6236 

21.7631 

70.0366 

- 43.4130 

2100 

6.3875 

4.9307 

26.9326 

22.0020 

67.0031 

- 40.0705 

2200 

6.4894 

4.9992 

27.2322 

22.2329 

64.2502 

- 37.0180 

2300 

6.5852 

5.0661 

27.5228 

22.4566 

61.7410 

- 34.2182 

2400 

6.6755 

5.1313 

27.8050 

22.6736 

59.4447 

- 31.6398 

2500 

6.7611 

5.1948 

28.0792 

22.8844 

57.3357 

- 29.2565 

2600 

6.8423 

5.2566 

28.3460 

j 23.0893 

55.3921 

- 27.0461 

2 700 

6.9198 

5.3168 

28.6057 

23.2889 

53.5954 

- 24.9897 

2800 

6.9939 

5.3754 

28.8587 

23.4833 

51.9297 

- 23.0710 

2900 

7.0650 

5.4324 

29.1054 

23.6729 

50.3814 

- 21.2761 

3000 

7.1336 

5.4880 

29.3460 

23.8580 

48.9387 

- 19.5927 

3100 

7.1999 

5.5422 

29.5810 

24.0389 

47.5912 

- 18.0102 

3200 

7.2642 

5.5950 

29.8106 

24.2157 

46.3300 

- 16.5194 

3300 

7.3267 

5.6465 

30.0351 

24.3886 

45.1472 

- 15.1120 

3400 

7.3877 

5.6968 

30.2548 

24.5579 

44.0357 

- 13.7809 

3 500 

7.4474 

5.7460 

30.4658 

24.7238 

42.9895 

- 12.5197 

3600 

7.5060 

5.7941 

30.6804 

24.8863 

42.0030 

- 11.3226 

3700 

7.5636 

5.8411 

30.8868 

25.0457 

41.0714 

- 10.1846 

3800 

7.6203 

5.8872 

31.0893 

25.2021 

40.1904 

- 9.1011 

3900 

7.6764 

5.9324 

31.2880 

25.3556 

39.3560 

- 8.0680 

4000 

7.7318 

5.9767 

31.4830 

25.5064 

38.5647 

- 7.0817 

4100 

7.7868 

6.0201 

31.6746 

25.6545 

37.8133 

- 6.1387 

4200 

7.8415 

6.0628 

31.8629 

25.8001 

37.0991 

- 5.2361 

4300 

7.8958 

6.1048 

32.0481 

25.9432 

36.4193 

- 4.3712 

4400 

7.9499 

6.1462 

32.2302 

26.0840 

35.7716 

- 3.5414 

4500 

8.0039 

6 . 1868 

32.4095 

26.2226 

35.1540 

- 2.7445 

4600 

8.0579 

6.2269 

32.5860 

26.3590 

34.5644 

- 1.9784 

4700 

8.1118 

6.2665 

32.7599 

26.4934 

34.0010 

- 1.2411 

4800 

8.1658 

6.3055 

32.9312 

26.6257 

33.4622 

- 0.5310 

4900 

8.2200 

6.3440 

33.1001 

26.7561 

32.9465 

0.1537 

5000 

8.2742 

6.3821 

33.2667 

26.8847 

32.4525 

0.8143 


5100 
5200 
5300 
5400 
5 500 

5600 
5700 
5 800 
5900 
6000 

6100 

6200 

6300 

6400 

6500 

6600 

6700 

6600 

6900 

7000 

7100 

7200 

7300 

7400 

7500 

7 600 
7700 
7800 
7900 
8000 

8100 

8200 

8300 

8400 

8 500 

8 600 
0700 
8800 
8900 
9000 

9100 

9200 

9300 

9400 

9500 

9600 

9700 

9 600 
9900 

10000 


TABLE m. - Concluded. THERMODYNAMIC PROPERTIES OF H 3 



c 

H„ - H n 


T 0 

R 

RT 

8.3287 

6.4157 

8.3835 

6.4569 

8.4385 

6.4938 

8.4938 

6.5303 

8.5495 

6.5665 

8.6056 

6.6024 

8.6621 

6.6381 

8*7191 

6.6735 

8.7765 

6.7086 

8.8344 

6.7436 

8.8929 

6.7783 

8.9518 

6.8129 

9.0114 

6.8473 

9.0715 

6.8816 

9.1323 

6.9158 

9.1936 

6.9458 

9.2556 

6.9838 

9.3183 

7.0176 

9.3816 

7.0514 

9.4456 

7.0852 

9.4601 

7.1185 

9.4733 

7.1511 

9.4849 

7.1830 

9.4948 

7.2142 

9.5030 

7.2447 

9.5096 

7.2744 

9.5146 

7.3035 

9.5180 

7.3319 

9.5198 

j 7.3596 

9.5201 

! 7.3866 

9.5189 

7.4129 

9.5162 

7.4386 

9.5120 

7.4636 

9.5064 

7.4879 

9.4993 

7.5116 

9.4908 

7.5347 

9.4809 

7.5571 

9.4697 

7.5789 

9.4571 

7.6001 

9.4432 

7.6207 

9.4280 

7.6406 

9.4114 

7.6599 

9.3937 

7.6787 

9.3746 

7.6968 

9.3544 

7.7144 

9.3329 

7.7313 

9.3102 

J 7.7477 

9.2864 

7.7636 

9.2614 

7.7788 

9.2352 

7.7935 


33.4311 

27.0114 

33.5934 

27.1365 

33.7536 

27.2598 

33.9119 

27.3815 

34.0682 

27.5017 

34.2228 

27.6203 

34.3756 

27.7375 

34.5267 

27.8533 

34.67 63 

27.9676 

34.8243 

28.0807 

34.9708 

28.1924 

35.1158 

28.3029 

35.2595 

28.4122 

35.4019 

28.5203 

35.5431 

| 28.6273 

35.6829 

28.7331 

35.8217 

28.8379 

35.9592 

28.9416 

36.0957 

29.0443 

36.2312 

29. 1460 

36.3653 

29.2467 

36.4977 

29.3465 

36.62 84 

29.4454 

36.7575 

29.5433 

36.8851 

29.6404 

37.0110 

29.7365 

37.1353 

29.8318 

37.2581 

29.9262 

37.3794 

30.0198 

37.4991 

30.1126 

37.6174 

| 30.2045 

37.7342 

| 30.2956 

37.8495 

30.3859 

37.9634 

30.4754 

3 8 . 07 58 

30.5642 

38.1869 

30.6522 

38.2966 

30.7394 

38.4048 

30.8259 

38.5118 

30.9117 

38.6174 

30.9967 

38.7216 

31.0810 

38.8246 

: 31.1646 

38.9262 

! 31.2476 

39.0266 

31.3298 

39.1257 

31.4113 

39.2235 

31.4922 

39.3201 

31.5724 

39.4155 

31.6519 

39.5097 

31.7308 

39.6026 

31.8091 


31.9789 

1.4522 

31.5246 

2.0687 

31.0885 

2.6651 

30.6696 

3.2422 

30.2669 

3.8013 

29.8796 

4.3432 

29.5069 

4.8687 

29.1480 

5.3788 

28.8022 

5.8741 

28.4689 

6.3553 

28.1475 

6.8232 

27.8374 

7.2784 

27.5381 

7.7214 

27.2491 

8.1528 

26.9699 

8.5731 

26.7001 

8.9828 

26.4393 

9.3824 

26.1871 

9.7722 

25.9430 

10.1527 

25.7069 

10.5243 

25.4780 

10.8873 

25.2556 

11.2421 

25.0395 

11.5889 

24.8294 

11.9282 

24.6250 

12.2601 

24.4260 

12.5849 

24.2323 

12.9030 

24.0437 

13.2144 

23.8598 

13.5195 

23.6806 

13.8185 

23.5057 

14.1116 

23.3352 

14.3990 

23. 1686 

14.6808 

23.0060 

14.9573 

22.8472 

15.2287 

22.6919 

15.4950 

22.5401 

15.7564 

22.3917 

16.0132 

22.2464 

16.2654 

22.1042 

16.5131 

21.9650 

16.7566 

21.8286 

16.9959 

21.6950 

17.2312 

21.5641 

17.4625 

21.4357 

17.6900 

21.3097 

17.9138 

21.1861 

18.1340 

21.0648 

18.3507 

20.9457 

18.5640 

20.8287 

18.7739 


TABLE IV. - COMPARISON OF CONSTANTS FOR THE GROUND ELECTRONIC STATE OF 


OBTAINED FROM THEORY AND EXPERIMENT 


Investigators 


Method 

Equilibrium 

Dissociation 

Fundamental 

Refer- 




internuclear 

energy, 

vibrational 

ence 




distance, 

D o> 

frequency, 





r e’ 

-1 

cm 

CD*, 

_1 





cm 


cm 


Cohen, Hiskes, and Riddell 

Theory 

(a) 

21343.71 

2190.96 

22 

Wind 


Theory 

1 . 057 X 10' 8 

21379.36 

2191.34 

15,23 

Hunter and Pritchard 


Theory 

1.061 

21379.22 

2191.23 

24 

Richardson 


Experiment 13 

1.057 

(c) 

2173 . 

45 

Beutler and Junger, 






f 46 

Herzberg and Monfils, 

r 

Experiment^ 

(a) 

21363 . 

(a) 

16 

and Namioka j 






Li« 

Kerwin, Marmet, and Clarke 

Experiment 6 

(a) 

(a) 

2194 . 

47 

Frost, McDowell, and Vroom 

Experiment^ 

1 < a > _.j 

(a) 

2178 . 

48 


a Not given. 

^Extrapolations from Hg bands. 
c No experimental value. 

^Obtained from dissociation and ionization energies of Hg. 
e Obtained from monoenergetic electron collision experiment. 
^Obtained by photoelectron spectroscopy. 
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TABLE V. - COMPLETE TABLE OF ROTATIONAL TERM VALUES OF GROUND ELECTRONIC STATE OF H+ 


Vibrational quantum number, v 


gular mo- 











mentum 

quantum 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 







1 




number, 

K 




Rotational term value, f, cm - 




0 

-0. 

-0. 

-0. 

0. 

0. 

-0. 

0. 

-0. 

-0. 

0. 

1 

58.26 

55.19 

52.21 

49.36 

46.56 

43.84 

41.10 

38.45 

35.79 

33.13 

2 

174.18 

165.11 

156.26 

147.73 

139.51 

131.01 

123.09 

115.01 

107.02 

98.95 

3 

347.14 

328.91 

311.25 

294.20 

277.61 

260.98 

244.99 

228.96 

213.05 

197.02 

4 

575.62 

545.29 

515.97 

487.63 

459.96 

432.55 

405.88 

379.31 

352.88 

326.27 

5 

858.04 

812.59 

768. 80 

726.53 

685.04 

644.28 

604.36 

564. 81 

525.36 

485.68 

6 

1191.99 

1128.80 

1067. 83 

1008.93 

951.15 

894.49 

838.91 

783.74 

72 8.74 

673.40 

7 

1575.00 

1491.61 

1410. 87 

1332.52 

1256.23 

1181.25 

1107.52 

1034.34 

961.32 

887.91 

8 

2C05.08 

1898.47 

1 795.45 

1695.72 

1598.02 

1502.42 

1408.17 

1314.65 

1221.24 

1127.15 

9 

2478.71 

2346.58 

2218.91 

2095.01 

1974.19 

1855.50 

1738.48 

1622.31 

1506.19 

1389.34 

10 

2592 .87 

2833.21 

2678. 57 

2528.60 

2382.13 

2238^21 

2056.29 

1955.26 

1814.18 

1671.78 

11 

3545.07 

3355.23 

3171.58 

2993.38 

2819.15 

2647.94 

2478.90 

2310. 01 

2142.67 

1972.92 

12 

4 131.31 

3909.69 

3 695. 04 

3486.52 

3282.60 

3082.04 

2683.9? 

2686.87 

2409.32 

2289.01 

13 

4 749.06 

4493.25 

4245. 87 

4005.26 

3769.73 

3537.95 

3308.99 

3080.56 

2 851.82 

2619.99 

14 

5 394.80 

5103.40 

4 821.22 

4546.76 

4277.70 

4012.95 

375C.92 

3489.67 

3227.28 

2961.42 

15 

6C65.55 

5736. 82 

5418.24 

5108.23 

4804.05 

4504.53 

4208.01 

3911.79 

3613.84 

3311.13 

16 

6758.13 

6390.63 

6034. 20 

5686.73 

5346.04 

5010.09 

4676.84 

4343. 74 

4008.29 

3666.66 

17 

7469.61 

7061. 87 

66 66. 09 

6280.13 

5901.04 

5527.05 

5155.59 

4783.69 

4408.36 

4025.39 

18 

8197.14 

7747.86 

7311.21 

6885.21 

6466 .65 

6052.98 

5641. 53 

5229.08 

4811.49 

4384.58 

19 

8537.96 

8445. 72 

7967.31 

7500.04 

7040.26 

6585.40 

6132.26 

5676. 96 

5215.31 

4741.35 

20 

9689.47 

9153.28 

8631.77 

8121.77 

7619.57 

7121.95 

6625.32 

6125.26 

5616.62 

5092.62 

21 

10448.80 

9867.93 

9302.22 

8748.36 

8202.41 

7660.38 

7118.57 

6571.27 

6012.95 

5434.81 

22 

11214.11 

10587.39 

9976.39 

9377.44 

8786.05 

8198.13 

7608.88 

7012.35 

6400.93 

5764.36 

23 

11582.93 

11309.45 

10652.09 

10006.03 

9368.73 

8733.10 

8094.45 

7445.72 

6777.61 

6075.87 

24 

12753.12 

12031.99 

11327.24 

10634.44 

9948.28 

9263.00 

8572.50 

7868.30 

7138.29 

6362.34 

25 

13522.70 

12753. 00 

11999. 80 

11258 .42 

10522.05 

9785.28 

9040.28 

8276. 75 

7478.30 

6610.26 

26 

14289. 72 

13470.59 

12667. 86 

11876.02 

11088.34 

10297.63 

9494.78 

8666.33 

7789.16 


27 

15G52.35 

141B2.92 

13329.47 

12485.95 

11644.71 

10797.32 

9932.47 

9032.08 

8054.29 


28 

15809.09 

14888. 16 

13982. 77 

13085.89 

12188.77 

11281.33 

10348. 72 

9364.56 



25 

16557.66 

15584. 59 

14625. 88 

13673.70 

12717.98 

11746.15 

10737.97 

9640. 79 



30 

17297.25 

16270. 72 

15256. 87 

14247.14 

13229.34 

12187.14 

11089.39 




3 1 

1 8C25.43 

16944.20 

15873. 75 

14803.69 

13719.23 

12597.55 





32 

18741.38 

176C3. 89 

1 64 74. 4 1 

15340.41 

14182.90 

12964.77 





33 

19443.30 

18247. 75 

17056. 51 

15853.73 

14612.85 






34 

2 0129.59 

18873. 78 

L 761 6. 98 

16338.61 

14993.43 






35 

2 C 798 . 63 

19479. 77 

18152.43 

16786.77 







36 

21448.63 

20062. 99 

18657.86 

17178.21 







37 

22077.51 

20620.03 

19124.20 ! 








38 

22682.94 

21145.94 









39 

23261 .79 

21631.45 









40 

23809.66 










41 

24319.10 





















Total an- 




Vibrational quantum number, 

V 




gular mo- 





















mentum 

10 

11 1 

12 

13 

14 

15 

16 

17 

18 

19 

quantum 

















-1 




number, | 
K 




Rotational tern 

value, f, cm 




0 

0. 

-0. 

0. 

-0. 

0. 

0. 

0. 

0. 

-o. 

0. 

1 

30.44 

27. 70 

24. 84 

21.93 

LB. 84 

15.54 

11.^2 

7.90 

3.42 

0.63 

2 

91 .06 

62.68 

74. 32 

65.48 

56.24 

46.27 

35.42 

23.28 

9.50 


3 

181.10 

164.62 

147.69 

130.06 

111.53 

91.59 

69.81 

45.37 

17.71 


4 

299.71 

272.33 

244.14 

214.77 

183.67 

150.53 

114.05 

72.90 

25.67 


5 

445.69 

404. 70 

362.48 

318.46 

271.78 

221.58 

166.74 

104.15 



6 

617.54 

560.32 

501.26 

439.66 

374. 15 

303.68 

226. 12 

136.67 



7 

813.62 

737. 50 

658. 87 

576.70 

489.36 

394.71 

289.95 

166.48 



8 

1C31.95 

934.43 

833.41 

727.72 

615.00 

492.44 

355.45 




9 

1270.75 

1149.06 

1022.89 

890.46 

748.99 

594.54 

418.72 




10 

1527.62 

1379.09 

1225.09 

1062.86 

888.81 

697.17 

474.36 




11 

1800.31 

1622.58 

1437.66 

1242.14 

1031.56 

797.23 





12 

2C86.50 

1876.79 

1658.01 

1425.91 

1173.91 

889.08 





13 

2383.67 

2139.26 

1883.20 

1610.65 

1311.80 

961.38 





14 

2689.31 

2407.35 

2110.86 

1792.90 

1439.82 






15 

3C00 .78 

2678.22 

233 7.2 5 

1968.69 

1548.68 






16 

3-15.39 

2948. 77 

2559.13 

2132.20 







17 

3630.21 

3215. 73 

2771.66 

2274.33 







18 

3542.16 

3475.33 

2969.24 








19 

4247.94 

3723.01 

3142.26 








20 

4543.56 

3952.84 









21 

4 824.31 

4154.04 









22 

5C83 .64 










23 

5310.09 
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TABLE VI. - INPUT FORM WITH DATA FOR SAMPLE PROBLEM FOR COMPUTER PROGRAM FOR WKB 


ROTATIONAL ENERGIES OF DIATOMIC MOLECULES AND MOLECULAR IONS 


[Top row contains card column numbers. ] 


10 


NL 


15 


13352.53 


JL ARRAY 
4660 6.72 
8904.35 


20121 


3013 ' 1 _ 


39310 .JL 2 53 107.5 2 J 27823.25 

5695 .57 [ 3^423_770 1 1866. 97__ 


40 141 


50 1 51 


60)61_ _ . 70J71 


80 


23313.34 

860.12 


19460.11 
278.37 _ 


16166. 

26.27 


64 


1 10111 20 [21 ‘ 30 P 51 40 141 

EL ARRAY ► 

.47625 " j .50271 j .52917 1 .55563 _ | .58208 

.71438 f -76729 | .82021 ] .87313 J_. 92604 


50 1.51 __ _ ‘ '60 1 6ll _ __70 1 71 80 

_r_.60854“_ j .6350.0 “ j .66146 JJ - 

I -97896 _ 1 1.03188 J 


1 10 

11 

NR 

UR ARRAY 

15 

30.24 

17689.41 

19994.73 


20 1 21 . 30 [31 ' ^40 [41 

j 391.0 4 “ L 153 4.71 ± 2934.51 

1 21249 .66 j 21891.94 J 22208.96 


50l51 " _60 1 61 

I 4753.17 ' 1 j 8100^62 
j 22363.11 j 22438.30 


70 [71 80 j 

’ - j 

_| 1381.8,30 1 

_ L22475. 68 J 


1 10)11. 

RR ARRAY 

1.08479 

1.16417 

3.20146 

5.67771 


20)21 __ 3013.1 40)4 1 „50[51 .. 6'0]jl * '1 . 70]71 ; l8jOJ 

T727 000 lV402 29~~ ' I i' .~55458~ j T. 77271 1~2. 24896 j 2.72521 ’ J 

1 4.15396 J 4.6 5021 . TT. 10646 . 1 5.58^7JL _J_ 6.05896 j J 


1 10 

ii 

20121 

30 [31 

40 [41 

_ 50 |51 _ 60J 61 

70 

NJ 

AJ ARRAY » 


5 

i. 

J 30. 

I 41.. 

. ±43 

1 44. 1 


1 5 |6 

NSTOP ARRAY 

72111 - 

„ _15 J 1 6 7 . 

“20l21 _25j2_6 

30 | 
1 

i 



1145 ) 'll 45 j " 1143 | 1143 [ _ . ' 1143J 


rr~. 

T2 

n 2 

NV 

GV ARRAY - 

18 

1145.08 

13932.54 

15299.56 

21679.25 

22095.95 


20J 21 


30 13 1 40141 _ "50 I5l_ 


60161 


70 1 71 _ 80 1 


3356.71 


16554. 96^ 
22371.29 


5400.56" 1 7540/29 _ 

_177 00.85 _187J3 . 55 _ 


9161.62 
_1 9650. 84. 


10861.50 

20450.59 


12450.33 
21128 . 60 


1 13 

14 26 

RE 

RM 

1.05685 

. 8364900E-24 


39J40 

T0L_ _ J N 


L_ 

lo 

t 

45 49 

50 

54 

55 ‘ 59 

60 64 

65 ' 69 

± NZ 

NMESH 


N 

NRITE 

N PUNCH 

M 

1 100 

li43 


10 

0 

0 

5[: 


• NAME 


75 
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TABLE VC. - OUTPUT FOR SAMPLE PROBLEM FOR COMPUTER PROGRAM FOR WKB ROTATIONAL 


ENERGIES OF DIATOMIC MOLECULES AND MOLECULAR IONS 


Descrip- 

tion 3 ^ 


C NL 15 UL 4 « 660872E 04 3.931012E 04 3.310752E 04 2.782325E 04 2.331334E 04 1.9460116 04 1.616664E 04 1.335253E 04 8.904350E 03 
5.695570E 03 3.423700E 03 1.866970E 03 8.601200E 02 2.783700E 02 2.627000E 01 

RL 4.762500E-01 5.C271Q0E-01 5.291700E-01 5.556300E-01 5.820800E-01 6.085400E-01 6.3500QCE-C1 6.614600E-01 7.143800E-01 
7.672900E-01 8.202100E-01 8.731300E-01 9.260400E-01 9.789600E-01 1.0318806 00 

NR 15 UR 3. C240C0E 01 3.94040CE 02 1-334 71 0E 03 2-934510E 03 4-753170E 03 8.1CC620E 03 1.381830E 04 1. 768941E 04 1.999473E 04 
2.124966E 04 2.189194E 04 2.220896E 04 2.236311E 04 2.243830E 04 2.247568E 04 

RR 1.08479CE 00 1.164170E 00 1.27000CE 00 1.402290E 00 1.534580E 00 1.772710E 00 2-24896CE 00 2-725210E 00 3.201460E 00 

3.677710E 00 4.15396CE 00 4.63021CE 00 5.106460E 00 5-582710E 00 6.058960E 00 

NJ 5 AJ l.COOOOOE CC 3.000000E 01 4.100000E 01 4.300000E 01 4.400000E 01 

NSTOP 1143 1143 1143 1143 1143 

NV 18 GV 1.145C80E 03 3.336710E 03 5.400560E 03 7.340290E 03 9.161620E 03 1.086150E 04 1-245033E 04 1.393254E 04 1.529956E 04 
1.655496E 04 1.770085E 04 1.873355E 04 1.965084E 04 2.045059E 04 2.112860E 04 2-167925E 04 2.209595E 04 2-237129E 04 


[re 

1.056850E 00 RM 

. 364900E— 25 TOL 

2.50Q000E-01 N2 

100 NHESH 1143 N 10 

NR IT E 

0 

NPUNCH 0 M 

5 

H2+ 

X 

''vv 

0. FV 58.02 

ERROR 

-0.24 

VV 1. FV 

55.17 

ERROR -0.00 


VV 

2. 

FV 

52.18 

ERROR 

_ 

0.25 

vv 

3. FV 49.38 

ERROR 

-0.08 

VV 4. FV 

46. 2*> 

ERROR -0.22 


VV 

5. 

FV 

43.73 

ERROR 

- 

0.01 

vv 

6. FV 41.12 

ERROR 

0. 

VV 7. FV 

38.50 

ERROR -0.00 


VV 

a. 

FV 

35.79 

ERROR 


0.00 

vv 

9. FV 33.05 

ERROR 

-0.17 

VV 10. FV 

30.50 

ERROR -0.00 


VV 

n. 

FV 

27.79 

ERROR 

- 

0.01 

vv 

12. FV 24. E5 

ERROR 

-0.00 

VV 13. FV 

2 1.98 

ERROR -0.10 


VV 

14. 

FV 

18. 77 

ERROR 


0.00 

vv 

15. FV 15.57 

ERROR 

-0.05 

VV 16. FV 

11.94 

ERROR -0.08 


VV 

17. 

FV 

7.94 

ERROR 

” 

0.17 

vm 17.65541 Fill 

22477. 

50 

V(IOO) -0.50000 FtlOO) 

60.08 









AJ 

1.00000 R 4 6 . 0589564E 

00 RRR (NSTOP) 6.0589564E 

00 5 

0 0 0 

0 


1 

0 

0 0 

0 

0 

0 

VV 

0. FV 17287. E5 

ERROR 

-0.00 

VV 1. FV 

16257.4 2 

ERROR 0.19 


VV 

2. 

FV 

15248.52 

ERROR 


0.01 

VV 

3. FV 14241.00 

ERROR 

0.01 

VV 4. FV 

13227.36 

ERROR 0.01 


VV 

5. 

FV 

12192.86 

ERROR 


0.01 

VV 

6. FV 11100.11 

ERROR 

-0.03 

VV 











V{ l ) 6.48722 Ftl) 

23699. 

72 

VUOO) -0.50000 F(IOO) 

17808.08 









AJ 

30.0000C R 4 4.26165326 

00 RRR l NS TOP) 6.0589564E 

00 5 

0 0 0 

0 


0 

0 

0 0 

0 

0 

0 

VV 

0. FV 243C8.23 

ERROR 

0.01 

VV 











Vll > 0.44822 Fll) 

25618. 

54 

V(100) -0.50000 FC100) 

25208.24 









Laj 

41.0000C R 4 3.2 549785E 

00 RRR (NSTOP) 6.0589564E 

00 5 

0 0 0 

0 


0 

0 

0 0 

0 

0 

0 


fv( 1 ) -0. 
\aj 43.0 


42914 F{ 1 ) 26218.21 VtlOO) *0.50000 F ( 100 1 26198.25 

00000 R 4 2.86 E4221E 00 RRR(NSTOP) 6.0589564E 00 5 0 0 


•s A J 44.0000C R 4 0. 


KRR(NSTOP) 6. 0589564E 00 


{ V( II -0.5000C Fill 
AJ 43.00000 R4 0. 


26300.12 VUOOI -0.50000 F<100) 26300.12 

RRR (NSTOP) 6. 0589564E 00 5 0 0 


0 0 0 
REO 00000 FIL« 


See p. 43. 
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TABLE Vm. - THERMODYNAMIC PROPERTIES OF Hg 


T, 

c p 

h t- h o 

9 ° 

S T 

<4-h 0 

h t 

g t 

°K 

R 

RT 

R 

RT 

RT 

RT 

298-15 

3.5224 

3.4623 

17.1098 

13.6475 

606.3907 

- 589.2810 

300 

3-5227 

3.4626 

17.1316 

13.6690 

602.6731 

- 585.5415 

400 

3.5509 

3.4806 

18.1482 

14.6676 

452.8884 

- 434.7402 

500 

3.6100 

3.5000 

18.9463 

15.4463 

363.0263 

- 344.0800 

600 

3.6947 

3.5251 

19.6117 

16.0865 

303.1304 

- 283.5187 

700 

3.7922 

3.5562 

20.1885 

16.6322 

260.3607 

- 240.1723 

800 

3.8918 

3.5920 

20.7013 

17.1094 

228.2959 

- 207.5945 

900 

3.9874 

3.6307 

21.1653 

17.5347 

203.3675 

- 182.2022 

1000 

4.0760 

3.6708 

21.5901 

17.9193 

183.4340 

- 161.8439 

1100 

4.1568 

3.7114 

21.9824 

18.2710 

167 . 1324 

- 145.1500 

1200 

4.2301 

3.7516 

22.3473 

18.5957 

153.5542 

- 131.2069 

1300 

4.2967 

3.7910 

22.6886 

18.8975 

142.0704 

- 119.3818 

1400 

4.3576 

3.8294 

23.00 93 

19.1799 

132.2316 

- 109.2223 

1500 

4.4136 

3.8665 

23.31 18 

19.4454 

123.7086 

- 100.3967 

1600 

4.4658 

3.9023 

23.5984 

19.6961 

116.2543 

- 92.6559 

1700 

4.5149 

3.9369 

23.8706 

19.9337 

109.6799 

- 85.8093 

1800 

4.5615 

3.9703 

24.1300 

20.1597 

103.8387 

- 79.7088 

1900 

4.6065 

4.0026 

24.3778 

20.3752 

98.6148 

- 74.2370 

2000 

4.6502 

4.0339 

24.6152 

20.5813 

93.9155 

- 69.3003 

2100 

4.6931 

4.0643 

24.8431 

20.7788 

89.6658 

- 64.8226 

2 200 

4.7355 

4.0938 

25.0624 

20.9686 

85.8044 

- 60.7419 

2300 

4.7775 

4.1227 

25.2739 

21.1512 

82.2805 

- 57.0067 

2 400 

4.8193 

4.1508 

25.4781 

21.3273 

79.0521 

- 53.5740 

2 500 

4.8607 

4.1784 

25.6757 

21.4973 

76.0836 

- 50.4080 

2600 

4.9017 

4.2054 

25.8671 

21.66 17 

73.3451 

- 47.4780 

2700 

4.9419 

4.2319 

26.0528 

21.8209 

70.8109 

- 44.7580 

2800 

4.9811 

4.2580 

26.2333 

21.9753 

68.4591 

- 42.2258 

2900 

5.0189 

4.2836 

26.4087 

22.1251 

66.2709 

- 39.8621 

3000 

5.0548 

4.3087 

26.5795 

22.2708 

64.2298 

- 37.6503 

3100 

5.0885 

4.3333 

26.7458 

22.4125 

62.3214 

- 35.5756 

3 200 

5. 1194 

4.3574 

26.9079 

22. 5504 

60.5334 

- 33.6255 

3300 

5.1473 

4.3809 

27.0658 

22. 6849 

58.8546 

- 31.7888 

3400 

5.1718 

4.4038 

27.2199 

22.8160 

57.2754 

- 30.0555 

3500 

5.1925 

4.4261 

27.3701 

22.9440 

55.7870 

- 28.4169 

3600 

5.2093 

4.4476 

27.5166 

23.0690 

54.3818 

- 26.8652 

3700 

5.2218 

4.4684 

27.6595 

23.1911 

53.0530 

- 25.3935 

3800 

5.2301 

4.4883 

27.7989 

23.3106 

51.7944 

- 23.9955 

3900 

5 . 234 ) 

4.5074 

27.9348 

23.4274 

50.6005 

- 22.6657 

4000 

5.2336 

4.5256 

28.0673 

23.5417 

49.4664 

- 21.3990 

4100 

5.2290 

4 . 54 2 8 

28.1965 

23.6537 

48.3875 

- 20.1910 

4 200 

5.2201 

4.5590 

28.3224 

23.7634 

47.3598 

- 19.0374 

4300 

5.2072 

4.5743 

28.4451 

23.8708 

46.3797 

- 17.9346 

4400 

5. 1905 

4. 5885 

28.5646 

23.9762 

45.4437 

- 16.8791 

4 500 

5.1702 

4.6016 

28.6811 

24.0794 

44.5490 

- 15.8679 

4600 

5.1465 

4.6138 

28.7944 

24.1807 

43.6927 

- 14.8983 

4700 

5.1198 

4. 6248 

28.9048 

24.2800 

42.8723 

- 13.9674 

4800 

5.0901 

4.6348 

29.0123 

24.3775 

42.0855 

- 13.0731 

4900 

5.0579 

4.6438 

29. 1170 

24.4732 

41.3301 

- 12.2132 

5000 

5.0235 

4.6517 

29.2188 

24.5671 

40.6043 

- 11.3856 



TABLE VIE. - Concluded. THERMODYNAMIC PROPERTIES OF Hg 


T, 

°K 

R 

h t - h o 

RT 

qO 

zz 

R 

g t- h o 

RT 

h t 

RT 

r° 

- 

RT 

5100 

4.9870 

4.6587 

29.3179 

24.6593 

39.9063 

-10.5884 

5200 

4.9487 

4.6646 

29.4144 

24.7498 

39.2344 

-9.8201 

5300 

4.9089 

4.6696 

29.5083 

24.8387 

38.5872 

-9.0789 

5400 

4.8679 

4.6736 

29.5996 

24.9260 

37.9631 

-8.3635 

5500 

4.8258 

4.6768 

29.6886 

25.0118 

37.3610 

-7.6724 

5600 

4.7829 

4.6791 

29.7752 

25.0961 

36.7796 

-7.0045 

5700 

4.7394 

4.6805 

29.8594 

25.1789 

36.2179 

-6.3585 

5800 

4.6954 

4.6812 

29.9415 

25.2603 

35.6748 

-5.7333 

5900 

4.6511 

4.6810 

30.0214 

25.3403 

35.1493 

-5.1280 

6000 

4.6068 

4.6802 

30.0992 

25.4190 

34.6407 

-4.5415 

6100 

4.5624 

4.6786 

30.1749 

25.4964 

34.1480 

-3.9730 

6200 

4.5182 

4.6764 

30.2488 

25.5724 

33.6704 

-3.4216 

6300 

4.4742 

4.6735 

30.3207 

25.6472 

33.2073 

-2.8866 

6400 

4.4306 

4.6700 

30.3908 

25.7208 

32.7580 

-2.3672 

6500 

4.3874 

4.6660 

30.4592 

25.7932 

32.3219 

-1.8627 

6600 

4.3447 

4.6615 

30.5258 

25.8644 

31.8983 

-1.3725 

6700 

4.3026 

4.6564 

30.5909 

25.9344 

31.4868 

-0.89 59 

6800 

4.2611 

4.6509 

30.6543 

26.0034 

31.0867 

-0.4324 

6900 

4.2204 

4.6450 

30.7162 

26.0712 

30.6976 

0.0186 

7000 

4.1803 

4.6386 

30.7766 

26.1380 

30.3191 

0.4576 

7 100 

4.1410 

4.6319 

30.8357 

26.2038 

29.9506 

0.8850 

7200 

4.1024 

4.6248 

30.8933 

26.2685 

29.5919 

1.3014 

7300 

4.0647 

4.6174 

30.9496 

26.3322 

29.2425 

1.7072 

7400 

4.0278 

4.6097 

31.0047 

26.3950 

28.9020 

2.1027 

7500 

3.9917 

4.6017 

31.0585 

26.4568 

28.5701 

2.4884 

7600 

3.9564 

4.5934 

31.1111 

26.5177 

28.2465 

2.8647 

7700 

3.9^20 

4.5849 

31.1626 

26.5777 

27.9308 

3.2319 

7800 

3.8883 

*.5762 

31.2130 

26.6368 

27.6228 

3.5903 

7900 

3.8556 

4.5673 

31.2624 

26.6951 

27.3221 

3.9402 

8000 

3.8236 

4.5582 

31.3107 

26.7525 

27.0286 

4.2021 

8100 

3.7925 

4.5489 

31.3580 

26.8090 

26.7419 

4.6161 

8 200 

3.7621 

4.5395 

31.4043 

26.8648 

26.4618 

4.942b 

8300 

3.7326 

4.5300 

31.4497 

26.9198 

26.1882 

5.2616 

8400 

3.7038 

4.5203 

31.4943 

26.9740 

25.9207 

5.5736 

8 500 

3.6758 

4.5105 

31.5379 

27.0274 

25.6591 

5.8788 

8600 

3.6466 

4.5007 

31.5808 

27.0801 

25.4034 

6.1774 

8700 

3.6221 

4.4907 

31.6228 

27.1321 

25. 1531 

6.4696 

8 800 

3.5963 

4.4807 

31.6640 

27.1833 

24.9083 

6.7557 

8908 

3.5)712 

4.4706 

31.7045 

27.2339 

24.6687 

7.0358 

9 000 

3.5469 

4.4605 

31.7443 

27.2838 

24.43*2 

7.3101 

9 100 

3.5231 

4.4503 

31.7834 | 

27.3330 ; 

24.2045 

7.5788 

9200 

3.5001 

4.4481 

31.8217 I 

27.3816 

23.9796 

7.8421 

9300 

3.4777 

4.4299 

31.8594 ; 

27.4296 

23.7593 

8.1002 

9400 

3.4559 

4.4196 

31.8965 j 

27.4769 

23.5434 

8.3531 

9500 

3.4347 

4.4094 

31.9330 | 

27.5236 

23.3318 

8.6012 

9600 

3.4141 

4.3991 

31.9688 

27.5697 

23.1245 

8.8444 

9700 

3.3940 

4.3889 

32.0041 

27.6152 

22.9211 

9.0830 

9800 

3.3746 

4.3786 

32.0388 

27.6602 

22.7218 

9.3170 

9900 

3.3556 

4.3684 

32.0730 

27.7046 

22.5263 

9.5467 

10000 

3.3372 

4.3582 

32.1066 

27.7485 

22.3345 

9.7722 
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Partition function, 





Internuclear distance (unequal side), cm 




2.5X10" 8 


2 - 


1.5- 



Internuclear distance (equal sides), R, Bohrs 


0 .5 1 1.5 2x10 " 8 

Internuclear distance (equal sides), cm 

Figure 2. - Potential energy of Hj for isosceles and equilateral triangle 
configurations (reproduced from fig. 13 of ref. 5). Inner ordinate gives 
internuclear distance R' in Bohrs (1 Bohr = 5. 29167xl0~ 9 cm) shown in 
the sketch. Inner abscissa corresponds to two equal internuclear dis- 
tances R (in Bohrs) also shown in sketch. Additional ordinate and 
abscissa scales in centimeters have been added. Potential energies were 
computed along lines A, J, K, M, N, Q, and S. Contours are lines of 
constant potential energy in Hartrees (1 Hartree = 219 475 cm" 1 ) and were 
faired in. To be consistent with figure 4 and rest of this report, 1.00 
Hartree must be added to all potential energies in this figure. 














Energy, hcU, hcU n , or hc^ 



Figure 6. - Typical potential energy hcU and effective potential energy hcU 1 ' for 
diatomic molecule or diatomic molecular ion (not to scale). For j = 0, hcU and hcU” 
coincide. Vibration -rotation energy eigenvalues hc^(4,20)and ho?(0, 0) are also 
shown, as well as hc*>{v^ ax , 20), which is not an eigenvalue. Four interpolation 
regions 1, 2, 3, and 4 are indicated just above abscissa. 



Figure 7. - Block diagram of program and subprograms for WKB rotational 
energies of diatomic molecules and molecular ions. Direction of call is indi- 
cated by arrowhead. Common storage is indicated by letter C. 
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